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1 .  INTRODUCTION 


The  feasibility  of  new  helicopter  designs  depends  to  a great  extent 
upon  an  analysis  of  their  flying  qualities.  Five  programs  have  been 
developed  to  carry  out  feasibility  analysis: 

These  programs  evaluate: 

1.  Helicopter  Stability  Derivatives 
including  static  stability  analysis 

2.  Solution  of  Longitudinal  Equations  of  Motion 

3.  Longitudinal  Motion  (Time  History) 

4.  Solution  of  Lateral  Equation  of  Motion 

5.  Lateral  Motion  (Time  History) 

Provisions  are  included  in  the  programs  for  investigation  of  stiMlit\ 
augmentation. 

A sample  calculation  is  included  for  a utility  helicopter.  The  he- 
licopter used  in  the  sample  calculation  was  chosen  to  demonstrate  the 
application  of  the  mathematics  only  and  the  type  of  evaluation  that  <.  an 
be  made  but  not  to  evaluate  present  capability. 

It  is  recognized  that  the  helicopter  chosen  is  not  a current  in- 
ventory item  but  was  used  because  the  dynamic  components  are  represen- 
tative of  present  design  and  there  was  extensive  flight  data  available 
with  and  without  the  SAS  to  compare  to  calculated  results. 

In  general,  the  procedures  and  notation  follow  those  established 


for  fixed  wing  aircraft.  Calculation  of  the  rotor  stability  derivatives 


are  added  and  the  familiar  aircraft  elevator,  aileron,  and  rudder  de- 
flections 5 and 

A’  R are  replaced  with  the  longitudinal  cyclic, 
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lateral  cyclic  and  tail  rotor  collective  Bj , , and  ©tR  Differences 

in  the  rotor  derivatives  between  articulated,  rigid  and  teetering  rotors 
mav  be  evaluated  with  the  stability  derivative  program. 


O 


2.  GENERAL  STABILITY  PROBLEM 

Normal  aircraft  have  six  equations  of  motion.  The  general  solution 
of  these  equations  is  a breakdown  of  two  sets  of  three  into  the  longitudinal 
and  lateral  modes.  Each  of  these  are  represented  by  a determinant,  the 
solution  of  which  is  a quartic  equation. 

There  are  four  possible  solutions  of  the  quartic  according  to  the 
roots  ,A,  being  real  or  complex  and  having  a positive  or  negative 
real  part. 

X = a i b i etc. 


A positive  real  part  indicates  an  unstable  or  diverging  motion  and 
the  complex  root  an  oscillation.  A measure  of  the  damping  of  its 
oscillation  is  the  time  to  damp  to  half  or  double  amplitude. 

T ._  .693 

l L = 

a 

ard  period 


2tt 


Another  common  way  to  write  the  stability  quartic  i s to  indicate 
the  natural  frequency,  u>  and  damping  ratio,  as 


The  characteristic  modes  ior  nearly  all  aircraft  in  the  longitudinal 
direction  are  two  oscillations,  one  the  long  period  or  "phugoid"  motion 
witi  very  light  damping  and  the  other  mode  a short  period  one  witn 
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vory  heavy  damping.  In  the  latexal  direction  there  are  three  normal 
modes.  Two  are  aperiodic,  one  a very  rapid  convergence  and  the  other 
a slow  divergence  known  as  "spiral  divergence."  The  third  mode  is  an 


* 


oscillatory  one,  the  well  known  "dutch  roll." 

In  general,  for  stable  oscillations,  it  is  desirable  to  have 
damping  ratios  on  the  order  of  0.6.  Specifications  for  various  types 
of  aircraft  generally  give  values  of  c,  ^ or  P and  t^  that  shall  be 
contractually  met,  for  example  see  MIL-H-8501A. 


S 


;— -I 


3 . STABILITY  EQUATIONS  REVIEW 


THE  MOMENTS  AND  FORCES  ACTING  ON 
A VEHICLE  ARE  WRITTEN  AS: 


AIRPLANE 

L-Cl1/2  AV2S 
D s CD 1/2  AV2S 
X = Cx  1/2  P V2S 
Y = Cy  1/2  P V2S 
Z = Cz1/2  P V2S 
L = Cl  1/2  P V2Sb 
M = Cm  1/2  p V2Sc 
N=Cn1/2  P V2 Sb 


HEUCOPTER 
= ClAVT|P2itR2 
s CpAVyip2  ir  R2 
s Cx  A Vtip2  it  R2 
s Cy  P Vyip2  it  R2 
s Cz  P Vjip2  ir  R2 
= Cl  P Vjip2  it  R3 
s Qa^  Vnp2  it  R3 
= CN^  VTIP2  w r3 


LIFT 

DRAG 

X AXIS  FORCE 
Y AXIS  FORCE 
7 AXIS  FORCE 
ROLLING  MOMENT 
PITCHING  MOMENT 
YAWING  MOMENT 
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The  equations  themselves  are  broken  down  into  two  parts,  the  left 
and  right  hand  sides.  The  right  hand  side  represents  the  inertia  forces. 


EF  -mu  EF  55  m u (£  + iji) 

x y 

EFZ  = m (w  - u S)  EL  = Ix  4> 


EM  = 


I 0 


EN  » 1 


z 


ip 


The  left  hand  forces  and  moments  are  series  expanded  in  terms  of 
the  changes  resulting  from  the  perturbations  in  the  linear  and  angular 
velocities  and  accelerations  of  the  aircraft.  These  are  known  as  the 
stability  derivatives  where 

3F  31-  . 3F 

E F = — ^ 6 + — <j>  + ^ I + etc. 

y 3 3 i ♦ 3 

It  is  the  purpose  of  this  note  to  develop  a program  for  estimating 
the  stability  derivatives  of  helicopters  and  additional  programs  to  utilize 
these  derivatives  in  solving  the  resulting  equationsof  motion.  A further 
objective  is  to  develop  step  by  step  solutions  so  that  a transient  motion 
of  the  helicopter  after  a disturbance  can  be  obtained. 

bach  of  the  forces  varies  with  the  disturbance  velocities  and  their  time 
derivatives: 


X- 

u 


u 


CX 


u 


DC 


x 


D a 


etc 
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The  standard  aircraft  longitudinal  equations  in  terns  of  the 
stability  derivatives  are: (for  example  see  reference  11). 


Xu  ♦ X a 
u a 


Z u ♦ Z a ♦ Z*a 
u a a 


Mu  + Ma  + M*  a 
u a a 


-we  + x*  6^  • m a 
6e  e 

z*e  * zt  & = n (w  - ve) 

e 6e  e 

M*e  ♦ m.  6 « i e 

e 6.  e y 


$e  » Bj  ■ Pitch  Cyclic  Longitudinal  Control  Deflection 


The  lateral  equations  are: 

v * V + ♦ 

L0e  ♦ i*u- 


Y.J+Y.  5n  * Yx 
* 6R  R {a  8 


H*  * + L6  6 a 

R a 


• • 

■ (v  + ViJO 


V 


V 


N;* 

y 


V 


♦ N 


5 6a  s V 
a 


* Tail  Rotor  Collective  Control  Deflection 

3 Lateral  Cyclic  Control  Deflection 
2 

/ q^R  Typical  Derivative  Coefficient 
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Solution  of  these  equations  takes  many  forms. 


For  static  lateral  stability  the  variation  of  control  deflection 
in  steady  sideslips  based  on  solution  of  the  determinants  become: 


-q  ttR  fC  fC  C 

dB  w 1 XB  1 6TR  nAl  ' "BTR  ~A1 


c c«  . 
1 J 


C0  rC  c C -C  . c rc  Cn  c„  C lw 

*6  t y0TR  nAl  ' 0TR  yAl  °B  ^ y0TR  V * ^OTR  yAl^y 


[ 6TR  A1  n0TR  A1J 


d9 


TR 


d<^> 


,C.  c C C,  dB 

t *B  nAl  ~ ng  *A1>  IT 
,C#  C C Cc  . 

C *0TR  nAl  ' n0TR 


dAl 

d<}> 


,C  C C C„  . dB 

( tg  n()TR  + nfi  eQTR)  ^ 

✓ Crt  L C Cn  V 

1 OTR  nAl  ' n0TR  A1J 


C Y 2 

Xg  = Vq»lT 


B n3 
qirR 


etc 


In  a steady  state  turn  B = 0 and  the  bank  and  turn  rates  become 

} . g -is*  e 


• t an<}> 
s g co?4> 


which  leads  to  the  result 
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TR 


•g  ^'.[e  (C1A1C‘;  - Cr.AICnj)  - cos,  (CSA1C*.  - CnA1  Cn;)] 
£ --  C C. 


\r  nAic- 


n\R  'Ale 


A1 


-«  V*  t • C'C  V % ♦ c W ' * C"eTi.C*ill 


C»  Cn 
0TB  "A1 
TR 


L 

cn  c* 
n0„n  A1 
TR 


In  pitch  the  variation  of  control  gradient  with  speed  is 


dBj 

air" 


C,  Cm  Cm  C, 

Zu  ma  • mu  Za 

^ _ 
zD1  m + m0,  z 
Blc  a Blc  a 


The  longitudinal  cyclic  pitch  control  required  in  a steady  turn  is 

g ISSUiSl  (c  C C C ) «•  f—  -1)  C,  C 
B.  » • v v z m;  - o z*'  lcos$  1 L m 
1 a 0 ot  0 r o a 


C2  \ - C»  CS. 

» Bj  c Bj 


or 


2 

E ,n  - 1 WC  C_  C-  C x . C.  C 

B,  . v )(  zo  *e  ' ",  *S>  * - » Lo  ", 


c c 

z„  an 
a B 


1 


W 


1 


C 2 

l0  ■ W/qirR  , n = Load  Factor 
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For  pitch  trim  in  stabilized  level  flight  the  load  factor,  n,is  1.  Moments 
about  the  C.G.  as  shown  in  Figure  1 give: 


X H CH  n2  A ___ 

R + R Cj  + 2Ct  p (2KfVf  - CL  TV)  - (2KVf  «r 


cLa  TV  I<t)2CT 


h 

R ' 2Ct 


(2KfVf 


For  steady  pull  ups  tho  pitch  damping  term,  Cm^,  may  be  added  as 
indicated  in  figure  1. 

Solution  of  the  static  stability  equilibrium  equations  has  been  added 
to  the  stability  determinant  program,  Appendix''^". 

Dynamic  stability  analysis  is  considered  in  separate  programs  as 
discussed  in  the  following  section. 


4.  LONGITUDINAL  ROOTS 

The  longitudinal  determinant  is 


C - — s 

X V 
II 


C7  + (C  - m')s 


a a 


C + C s 
m m* 

a a 


(C  + m')s 
C s - Iy's2 


where 


m'  = Sr  » l'y  = Iyy/qADROTR 
4 b 
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LONGITUDINAL 


4 3 2 

This  gives  the  stability  quartic  A's  + B's  + C's  + D's  + £'■ 


Let  A = 
Then: 

C , B = C_  , C 
x * Z 

u u 

= C , E = 
x ’ 
a 

cz  , D = cz. 

a a 

m' , F = Cz<  + m' 

a 

A' 

B' 

C' 

D' 

E' 

~ Diy 

-ADIy ' 

-AEIy ' 

AECm* 

9 

-BC.  C 
L m 

a 

m ' r , , 

rEly 

+ADCm* 

0 

-AFC 

m 

a 

EC,  C 
L mu 

v m* 

0 

-AFCm* 

a 

-BCC»6 

m'  cr 

+ — FC 
v m* 

a 

- - EC 
v m* 
6 

-BC  C 
L m* 
a 

+ -1  FC 
v m 

a 

CFG 

mu 

+BCIy ' 

DC.  C 
L mu 

These  are  solved  in  the 

program  of  Appendix  "B". 

For  hover  the  Z 

• • 

equation  Vy  term  becomes  ui. 

This  is  accomplished  by  setting  v = 1 in  the 

input  of  the 

program. 

It  is  easy  to  introduce  various  autopilot  gains  into  the  quartic.  For 
example, C52  = B^e  adds  these  terms  to  the  quartic  coefficients  (Table  B.l). 


AC' 

ad* 

AE' 

•f  DC52CR‘B1 

ADC52  CmBl 

AEC52CmBl 

m 1 

i-C.-Citl'C, 

V 52  a 2Bi 

-AC52Cm;CZB1 

-AC  .,Cm  C_ 
a L 

-m 1 

2i-EC.„CmD1 
v 52  B1 

~BCC52CmBl 

f 16 

! 
t 
i 

i 


AC' 


AD* 


AE' 


f - C52  Cm«CZ 


B1 


+CC  Cc,  C, 
mu  52  Z 


B1 


*C«  Cm*Cc~C 
Z ot  52  xD1 

U 


-C_  Cm  Cc_C 
Z a 52  xD1 
y B1 


+Cm  D C..C 
y 52  X 


B1 


♦C  C_  C._C 
my  Z 52  xD1 
a B1 


If  for  a given  aircraft  and  flight  condition  a set  of  desired 

longitudinal  phugoid  - short  period  characteristics  are  required,  a 
desired  quartic  is  known: 


A's4 

+ B's3  ♦ C's2  +D's  + E' 

Desired 

A'.s4 

+ BJ  s3  + C|s2  + D|s  + E' 

Aircraft  Base 

AA^s4  + AB's3  + AC's2  + AD's  + AE'  Increment 


For  each  flight  condition  it  is  possible  to  solve  for  a set  of 
feedback  gains  to  satisfy  this  condition.  Five  sets  would  be  required. 
For  example  with: 

c52  * B1/e  c*15  - Bi/ct 

C*51  = Bl/e  C'16  = Bl/a 

C'50  = Bl/0  C'17-B1/o 


From  the  expansion  of  the  stability  quartic,  obtain  the  coefficient 

increments  due  to  each,  see  Table  B.l. 

= AA  due  to  C'15  etc. 
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which  gives  for  each  flight  condition  the  set: 

A15‘C' 15+A16-C' 16+A50*C' 50+A51 •C'51+A52>C '52  = AA_ 

K 

B15-C'15+B16-C'16+B50.C'50+B51*C‘51+B52-C,S2  = ABR 

C15*C' 15+C16-C' 16+C50<C*  50+C51 •C'51+C52*C'52  = aCr 

D15«C' 15+D16’C' 16+DSO‘C' 50+D51 •C'51+D52*C'52  = ADR 

D15*C' 15+E16»C' 16+E50*C'50+E51 •E'51+E52*C'52  = A0R 

Solution  of  this  set  will  give  the  desired  stability. 

Scheduling  of  5 multiple  gains  through  the  entire  flight  regime 
becomes  quite  cumbersome.  Usually,  a hunt  and  trial  method  using  only 
one  or  two  gains  to  obtain  an  approximate  solution  is  used.  To  avoid  this 
problem  the  idea  of  the  high  gain-model  following  autopilot  has  been  used. 


Cgf%Z 


— J2}« 


'5* 


In  this  figure,  if  the  autopilot  gains  are  adjusted  so  that  a desired 
damping  ratio  of  .5  or  .6  is  obtained, and  if  the  variable  gain  is  made 
large  enough,  the  airframe  dynamics  will  tend  to  be  masked  and  the 
damping  ratio  set  in  to  the  autopilot  will  be  obtained. 

This  arrangement  is  similar  to  the  adaptive  autopilot  scheme 
(references  6,  10,  and  12), except  that  the  aircraft  function  is  multiplied 
by  the  autopilot  function  in  the  equivalent  circuit. 


18 


With  this  simplified  arrangement  an  insight  into  the  dynamic  stability 
problem  can  be  obtained  by  using  the  standard  quartic  root  procedure  as 
used  for  the  airframe  alone,  see  Table  B.l.The  required  inputs  have  been 
set  into  Appendix  B.  A corresponding  transient  analysis  is  set  up  with 

provision  for  a variable  gain  in  Appendix  C. 

•# 

Inspection  of  the  0 equation  shows  in  another  way  the  influence  of 
the  autopilot  terms 


6 ■ 


C u+C  a+C  a ♦ GVC52C  6 + (GVC51C  + C )0 
» m m*  mD1  m ' 

U a a B1  B1  0 


Iy'  - GV  * C50  C 


m( 


B1 


Thus  it  can  be  seen  in  the  denominator  an  increase  in  GV  or  C50 
is  equivalent  to  a change  in  the  moment  of  inertia. 
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5.  Lateral  Roots 


The  lateral  determinant  is 


where 


mV 


*xx  T , _ XYY 
> 


m'  ■ j , I*  - ? » \ 

qllR2  X qlir  Y qI!R 


-m's  + 


(Cy.  - m' )s 


ci*s 


-i;  s2  * c -s 

Z n>J> 


CL  + (CY*  + m'a)S 
-!xs2  ♦ C£-s 


+C  *s 

n<t> 


which  gives  the 


stability  quartic  A's^  + B*s3  + C's2  +D's  + E* 


A' 

B' 

C' 

D* 

E’ 

*'vxlz 

-C  1*1* 

Vlz‘x 

p 

CY  lZCl' 
B ♦ 

‘CttCYiCni 

C n c c. 

h % L 

-,,z % 

Cv  C .11 
YR  n^  X 

%"'Cn* 

-c  Co  c. 
ne  % L 

-'h% 

'm'C^Cn|  'C£BIZCL 

ra’Cnft 

C£BCn^CY4> 

~C€eIZCY4> 

SnBCY^. 

_Cn  BTXCY^ 

-C  . m 1 C p 
ng  <~1 
<P 

Cn?'h 

•-C  .Cp*  Cv» 
ng  fy  Yij> 

W4 

"cybc4c4 
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These  determinants  are  solved  in  the  program  of  Appendix  D.  For 
hover,  the  side  force  equation  term  Vifi  disappears.  This  is  accomplished 
by  setting  V«1  in  the  input  of  Appendix  D.  As  for  the  longitudinal 
case,  a simplified  autopilot  study  is  included.  Autopilot  gain  coef- 
ficients are  shown  in  Table  D.l. 

The  corresponding  transient  analysis  is  included  in  the  calculation 
of  Appendix  E. 
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6.  STABILITY  DERIVATIVES 

Helicopter  stability  derivatives  are  influenced  by  the  type  of 
rotor.  In  general  three  types  of  rotor  are  commonly  used. 

1.  articulated  rotor 

2.  rigid  rotor 

3.  teetering  rotor 

The  first  two  transfer  cyclic  control  moments  directly  to  the 
rotor  shaft.  The  cyclic  lift  couple  for  the  first  has  an  arm  equal  to 
the  radius  of  the  flapping  hinge  and  the  second  equal  to  the  radius 
of  the  blade  spanwise  center  of  lift.  The  teetering  rotor  develops 
its  moment  by  tilting  the  lift  vector  about  the  C.G.  No  moment  is 
developed  at  zero  lift. 

Two  programs  have  been  developed  for  estimating  stability  deriv- 
atives. The  first  program  was  based  on  combining  the  equations  and 
data  of  reference  1 and  2.  In  general,  this  work  was  based  on  the 
articulated  rotor  only.  In  addition  contributions  of  the  rotor  to  some 
of  the  small  damping  turns  were  neglected. 

For  this  reason  a second  program  based  on  a rotor  strip  analysis 
was  developed  to  calculate  complete  sets  of  stability  derivatives 
including  a complete  model  build  up.  The  first  program  is  shown  in 
table  A1  and  the  second  table  A2  of  appendix  A.  The  method  of  A1 
being  a closed  solution  is  much  cheaper  to  operate  than  program  A2 
which  requires  a great  deal  of  machine  iteration  and  thus  machine  time. 

In  general  values  obtained  from  these  methods  may  be  used  where 
wind  tunnel  o<-  other  test  data  are  not  available.  Some  useage  of 


these  methods  and  results  obtained  are  presented  in  the  following 
examples. 

7.  EXAMPLE  CALCULATION 

Several  initial  studies  have  been  completed  using  the  methods 
of  this  report.  These  are  shown  in  appendix  F. 

The  first  study  was  completed  using  the  abbreviated  derivative 
calculation  of  table  Al.  The  purpose  of  this  study  was  to  illustrate 
the  influence  of  variations  of  the  stability  derivatives  in  the  dynamic 
stability  of  a marginally  stable  basic  helicopter  examples  I and  II. 

In  addition,  the  static  stability  relations  of  a utility  helicopter 
was  investigated  in  example  III. 

Stability  derivatives  are  shown  in  Table  F.l,  Appendix  F for  typical 
study  helicopters  whose  characteristics  are: 


w 

Weight 

lbs. 

I 

11867 

II 

11867 

III 

14800 

CDA 

Helicopter  drag  area 

sq.  ft. 

36.5 

36.5 

21 

V 

Velocity 

ft/sec 

100 

0 

169 

23 


I 

II 

III 

TIPMS 

Main  Rotor  Tip  !!ach  No. 

.6 

.6 

.65 

SIGMA 

Main  Rotor  Solidity 

B*C/ 

(3. 14*R) 

.062 

.062 

.09 

ROTR 

Main  Rotor  Radius 

ft. 

28 

28 

28 

Cl 

Speed  of  Sound 

ft/sec 

1115 

1115 

1115 

RHO 

Density 

slugs/ 
ft  3 

.0024 

.0024 

.0024 

XR 

Distance  Rotor  to  CG  Horiz. 

ft-CG  fwd 

1 -.1 

-.1 

.5 

HR 

Distance  Rotor  to  CG  Vert. 

ft+CG 

below 

8.0 

8.0 

6.25 

HF 

Munks  KF  for  Fuselages 

.8 

.8 

.7 

VOLF 

Munks  Fuselage  Volume 

cu.  ft. 

1000 

1000 

1500 

ST 

Area  Horiz.  Tail 

sq.  ft. 

12.4 

12.4 

30 

SVT 

Area  Vert.  Tail 

sq.  ft. 

20 

20 

25 

LT 

Tail  Length 

ft. 

28 

28 

33.3 

TI 

Horiz.  Tail  Incidence 

rad. 

0 

0 

0 

SIGMA T 

Tail  Rotor  Solidity 

.167 

.167 

.167 

TRR 

Tail  Rotor  Radius 

ft. 

4.7 

4.7 

5.0 

DS 

Angle  Shaft  to  Body  Axis  (+<5  AHD  Z) 

Rad. 

.174 

.174 

.087 

SG 

1 for  DP=RHO*VTIP*VTIP , 
0 for  DP-RH0*V*V/2 

0 

1 

1.0 

BMF 

Blade  Mass  Factor  ■ C*RHO*5.73 
(R0T0R**2)/I 

10 

10 

10 

ETA 

Load  Factor 

1.0 

1.0 

1.0 

From  the  stability  derivatives  of  Table  F.l,  dynamic  flight 
characteristics  have  been  calculated  for  examples  I and  II.  These 
are  shown  in  Appendix  F,  Tables  F.2  through  F.4.  The  effects  of  changes 
in  stability  derivatives  and  autopilot  gains  are  shown.  Figures  F.l  to 
F.3  show  typical  motions  of  the  aircraft  with  and  without  typical  auto- 
pilot data. 

Finally,  static  stability  data  are  tabulated  in  Table  F.4  for 
example  III,  The  data  include  pitch,  equilibrium  sideslip  trim  and 
steady  turn  data.  These  data  are  also  plotted  in  Figures  F.4  to  F.8. 
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Figure  F.6  indicates  that  the  utility  helicopter  of  example  III  should 
have  adequate  longitudinal  control  for  a CG  range  of  16.0  inches.  This 
range  could  be  extended  to  27.0  inches  by  coupling  the  horizontal  tail 
to  the  longitudinal  cyclic  control.  However,  this  would  be  the  maximum 
amount  possible  since  the  gradient  of  control  position  'versus  speeds 
becomes  critical  at  the  aft  CG,  see  Figure  F.5. 

Some  recent  flight  data  is  available  on  a series  of  Bell  Hellicopters 
with  stability  augmentation  off.  These  helicopters  use  a 44  foot  diam- 
eter rotor  with  27  inch  chord  operating  at  tip  speed  of  745  feet  per 
second.  Aircraft  using  this  rotor  include  the  UH-1C  and  AH-1  series. 

Tests  of  the  AH-1  series  found  the  airplane  with  stability  augmentation 
off  to  be  dynamically  unstable  although  not  uncontrollable  for  several 
flight  conditions. 

Stability  derivatives  were  estimated  for  these  aircraft  at  hover 
at  7100  lbs  and  level  flight  at  40  and  140  knots  forward  speed.  All 
derivatives  were  estimated  by  the  rotor  strip  analysis  of  table  A-2. 

In  addition  for  the  hover  condition  derivatives  were  estimated  by  the 
method  of  both  tables  A1  and  A2.  These  are  tabulated  in  Table  F.6.  The 
estimated  differences  between  using  an  equivalent  teetering,  articulated 
or  rigid  rotor  at  hover  are  also  shown  in  Table  F.6. 

Calculated  values  of  the  period,  P,  and  time  to  damp,  t,  using 
these  derivatives  are  compared  with  flight  values  as  follows: 
where  P/t  +t  time  to  double 

-t  time  to  half 
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LONGITUDINAL  (NO  SAS) 


Calculated 

Flight 

AH-1 

Hover 

16/ Divergent 

Model 

40  kts 

5. 5/ +63 

16/+4 

(R) 

S/Neutral 

(Q) 

140  kts 

DEAD  BEAT 

DEAD  BEAT 

(Q) 

LATERAL  (NO  SAS) 

Hover 

5. /Divergent 

40  kts 

Divergent 

Divergent 

(R) 

5./-S 

(Q) 

140  kts 

2.8/100 

4.5/Neutral 

(Q) 

In  summary  the  methods  of  computing  stability  derivatives  shown 
in  table  A-l  and  A-2  both  need  continual  correlation  and  updating. 

The  methods  may  be  used  as  a basis  of  correlation  or  for  use  in  prelim- 
inary estimates  where  wind  tunnels  or  other  test  data  is  not  available. 
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APPENDIX  A 


HELICOPTER  STABILITY  DERIVATIVE  PROGRAM 


ooooooonnoonoononoonoonoooooorioriooooono 


TABLE  A-l 


CLOSED  SOLUTION  ANALYSIS 
HELICOPTER  STABILITY  DERIVATIVES 
C UTILIZATION  PROGRESS  REPORT 


THIS  PROGRAM  EVALUATES  STABILITY  DERIVATIVES  FOR  A SINGLE  ROTOR 
HELICOPTER.  THIS  FIRST  CUT  PROGRAM  IS  LARGELY  BASED  ON  METHODS 
OUTLINEO  IN  SECKELS  HOOK  COMBINED  WITH  DATA  FROM  SIKORSKY  AND 
BOEING  RC TOR  CHARTS. 

THE  REQUIRED  INPUT  CARDS  AND  PROPER  COLUMNS  ALONG  WITH  THE  SYMBOLS 
AND  IDENTIFICATION  ARE  PRESENTED  IN  THE  FOLLOWING  TABLE. 


CARO 

I 

COLUMN 

WO 

SYMB  L 
W 

WEIGHT 

LBS 

l 

ll-?C 

CDA 

HELICOPTER  DRAG  AREA 

SO.  FT. 

1 

21-3C 

V 

VELOCITY 

FT/SEC 

1 

31-AO 

TIPMS 

MAIN  ROTOR  TIP  MACH  NO 

1 

61-50 

SIGMA 

MAIN  ROTOR  SOLIDITY 

B*C/ ( 3. 16*R ) 

l 

51-6C 

ROTR 

MAIN  ROTOR  RADIUS 

FT. 

I 

61-  7C 

Cl 

SPEED  OF  SOUND 

FT. /SEC. 

2 

l- 10 

RHO 

DENSITY 

SLUGS/FT  3 

2 

11-20 

XR 

DISTANCE  ROTOR  TO  C.G.  HORIZ 

FT  -C.G.  FWD 

2 

21-3C 

HR 

01  STANCE  ROTOR  TO  C.G.  VERT . 

FT  ♦C.G.  BCL Ow 

2 

31-AC 

HF 

MUMKS  KF  FOR  FUSELAGES 

2 

61-SC 

VOLF 

MUNKS  FUSELAGE  V)LU*E 

CU.FT. 

2 

51-6C 

ST 

AREA  HORIZ  TAIL 

SQ.FT. 

2 

61-70 

SVT 

AREA  VERT  TAIL 

SQ.FT. 

3 

WO 

LT 

TAIL  LENGTH 

FT. 

3 

H-2C 

T I 

HORIZ  TAIL  INCIDENCE 

RAD. 

3 

21-  30 

SIGMA 

TAIL  ROTOR  SOLIDITY 

3 

31-AO 

TRR 

TAIL  ROTOR  RADIUS 

FT. 

3 

61-5C 

OS 

ANGLE  SHAFT  TO  BODY  AXIS  *S  AHD  Z 

RAD 

A 

1 — 1 0 

5G 

l FOR  OP  *RHO*VTI P*VT I P 0 FOR  DP 

=RH0*V*V/2 

A 

1 1-20 

HMF 

BLADE  MASS  FACT0R=C*RH3*5.73*(R0T0R**6)/ l 

A 

21-3C 

FTA 

LOAD  FACTOR 

A 

31-AO 

b r i 

RATIO  DBI/DTI  PITCH  CYCLIC  TO  TAIL 

INC  I QENCF 

11867. 

36.5 

TYPICAL  INPUT  DATA  ARE 
IOC'.  .6  .062?  28. 

1115. 

.0026 

-1. 

8. 

.8  1000.  12.6 

20. 

28. 

0. 

.167 

6.67  .176 

1. 

10. 

1. 

0. 

CALL  TANNRHUIPMS, SIGMA, CTOSH.CUOSH.CTSTH.O) 

CALL  TANNESI TIPMS,OMUHS,CLOS ,CDOSS ,CQOSS,CQOOSS , 0, ACSS ) 

77  REA0<5,m,CPA,V,  TIPMS,SIGMAfRCTR,Cl 
READ! 5,  lHH.,XR,HR,HF,VOLF,ST,SVT 
RE  ADI S, l ) TL, TI ,SIGMAR,TRR,DS 

C TO  CALCULATE.  LONGITUDINAL  CYCLIC,  BIC,  IN  PULL  UP  SET  6T A 

C (LOAD  FAC  TCR  ) IN  CARD  6.  NORMALLY  ETA  =1.  FOR  HOVFR  ETA  =1. 

REAOl 5, 1 I SO.^MF,-  TA.BTI 
WRITF (6,M  ) 

A- 2 


nnnno 


WRITE(6,2) 

2 FORMAT; iuXtl OH INPUT  OATA) 

WRITE!6,1>W,C0A,V,T!PMS,SIGMA,R)TR,CI 
WR ITE ( 6, 1 |RH0»XR«HR  *HF » VOLF , ST  »SVT 
WRITE (6, 1)TL,TI,SIGMAR,TRR,DS 
WRITE (6,  l)SG,RMF,ST4,BTI 
IF!ETA-1.)13»13»12 
12  ETB*22.2*1ETA-1.)/V 
GOTOIA 
li  ET8=0. 

14  CONTINUE 
51  FORMAT! IH1 ) 

VTIP*TIPM$*CI 

0P«RH0*V*V/2. 

OMUHS«V/VTIP 
AS4=SIGMA*5.  13/4. 

AD«3.i41*R0TR*R0TR 

CT»W/!AO*VTIP*VTIP*RHO) 

ELQ=-l./!!2./CT)**.5) 

CLOS*CT/SIGMA 
CT  AS  = 2.*CT/ ! 5. 73*S IGMAI 
CDQ$  = (C04*f)P )£CT„S/W 
IF ( OMUHS ) 60, 61 » 60 

60  CALL  TANNEStTIPMS, OMUHS, CLOS,COOS, COOS, CQOOS ♦ 1 * ACS2 ) 

GOT''  62 

61  CALL  TANNRH(TIPMH, SIGMA, GLOS  ,CQOS  ,CTSTH,l) 

OP»RHO*VTIP*VTIP 

DP1*DP 
WR ITE ( 6* 15 1 

62  CONTINUE 

XU*!OX/f>U)/M  * ( G/C T ) * ( fiR ) SECKEL  SYSTEM 

CXU  NATURAL  AIRCRAFT  SYSTEM* ( M/Q*AD » *XU  ALPHA*W/V 

CXU  NATURAL  HELICOPTER  SYSTEM  • USE  QM*RHO*VTIP*VTIP  V*l.  THEN  AL*W 

THEREFORE  SET  SG=l  FOR  HELICOPTER  SYSTEM,  SG*0  FOR  AIRCRAFT  SYSTEM 

NOT?  COR  CT  PM  IN  DYNAMIC  PROGRAMS  TO  MATCH 

IF ( SG ) 91 »Sl ,90 

90  DP=RHO*VTIP*VTIP 
WR I TE ( 6, 93 ) 

G0T092 

91  WR ITE! 6, 94) 

92  CONTINUE 
DPF=RH0*V*V/!2.*DP) 

CQA5>s2.*CQ0S/5.73 

Ql  = ! 4.2*  38 .*CTAS*?.* OMUHS )/ 1 0000. 

Q3»(-1.8-2.86*0MUHS**2. 51/10000. 

Q2S! CQAS-Qi )/!CTAS**2) 

ELS=-3.*  MUHS-! .92-.5*0MUHS)*02 
Q2=Q2+Q3*EL  S/CTAS 
ELW=£LS*CTAS 
TSRP=0MUHS*!2./CT)**.5 
IF(TSRP-l.)2C,20,2l 
?C  0L0W=.5*.25*TSRP 

OLOCT=( .25*.03*TSRP)/EL0 
OLOU= . 35*  TSRP 
Ell=.4*TSRP 
GOTO  28 

21  IF1TSRP-?. >22,22,23 

22  DLCW=.75+.2*(TS9P-1.) 

OLOC.T  = ( .2P-t  .06*!  TSRP-l.  ) ) )/ElO 
OLDU=. 35-1.14*! TSRP-l.M 


El  l=.4-<  .OR*(TSRP-l.  ) ) 

GOTO  2B 

23  IFITSRP-  3.)24t24,25 

24  0L0W*.95+(.25*( TSRP-2. ) » 

DLOCT=( .22-1 .04* ( TSRP-2. ) H/ELO 
DL0U=.2l- ( « l* ( TSRP-2  • )) 

ELl=.32-(  .04MTSRP-2.H 

GOTO  28 

?5  IFITSRP-6. >26,26,27 

26  OL0W=l. 

DLOCT*!.  18-U03  3 1MTSRP-3. > ) ) /FLO 
DLDU=.ll-( .03  33* < TSRP-3. ) ) 

EL  1= .28- ( , 04* ( TSRP-3 • ) ) 

GOTO  2 R 

27  DLCW=l. 
niOCT=.OH/Fi: 

DLDU=.Ol 

EL  1= . 16 

28  DL OUsOLOU+ . CC  5* AC  S2 

DCTOT  = U.  + l.ti*0MiiHS<'0HtJHS)/»3.*(l.-ASA*nLnCTn 
0CT0W  = .5*0L0Vi/(  l.-AS4*0LDCT) 

DCT0U=<3.*0MUHS/< l.+l.5*0MUHS*0MUHS) )*tCTAS~6LW/2.l 

OCTfHJ*!  DCTOU+  . 5*DL0U ) / ( 1 .-AS4*DLBCn 

R-ROTR 

31  UP=< ( XR/K ) ♦ ( 0MUMS*2 .*ELS*IHF*VOLF-l.43*TL*ST)/(5.73*SIGMA*AO*ROTR) 
1-0MUHS*CMUHS*  (HF* V0LF*0. *1 . 43*TL*ST*TI  )/ICT*A0-rtOTR))  ) 

0N  = (HR/R»-0MUHS*^MUHS*(HF*V0LF-l.A3*(BTI*l.  »*ST*TU/(CT*AO*R'.TR) 

H4  = I6.9'i*CMUHS**.92)-l.l*0MUHS*EL$ 

UP»UP*(  .v  l25*0MUHS/( U.4*CTASl*H4*CTAS-16.*ETe*R*(l.*ELS/4. )/(BMF* 
1VTIP ) )*H*/R 

bic=up/0-; 

TlC*(CTAS-FLW/2.  »*3.**>7.3/U.*l  ,5*0PUHS*OMUHS> 

H4P=1.*< IC.*<  l.-0MUHS)**2.9*>>*ELS 
H1P=1«S.7*CMUHS-.6R5*?L$*OMUMS**.60S 
H2P  = -.9*0MUHS*.0C  3*ELS 
H3P=8.-5.*0MUHS-( 1.5-2 . 25*0MUHS ) *EL  * 

H5P=l.0*3.l*GMUHS**2.04M ,25-.07*0MUHS)*ELS 
DCH0T  = (-  BIC*CTAS*(H1P+-2.*AS4*H2P*DL0CT))*DCT0T 
0CH0W*(-8lC+CTAS*(HlP+2.*H2Pn*f)CT0W 

DCH0U-(-8iC*CTAS*(HlP*2. * AS4*H2P*DLDCT ) ) *DCT  U*H2  P*CT  AS^DLIHJ  *H3P*CT 
1AS*CTAS*.C125/11,.46 
CCHf)Q*-l 6 ,*C  f AS*H5P/BMF 
nCHOBl=-  MUHS*Of.HOW-CTAS 
DCTnfil=-DCTOW*0MUHS 
DCTnP={0MUH$/4.)/(i,-AS4*0LDCT) 

AS5*AS4*SIGMAR /SIGMA 
C01=W/{VT!P*CTAS*DP*A0) 

CXUR=C01*(-nCH0U+0S*(DCH0W+0CTDUl-DS*DS*0CTDW) 

CXU=CXUR-2.*RH0*V*CDA/(DP*AD) 

CXW=COl*I-PCHDW+DS*OCTDW) 

CXA=CXW*V 

C02=W/(CTAS*DP*AD) 

CXTC=CO?*(OCTOT*DS-nCHOTJ 

CXBl=C02M0CTnBI*DS-DCHnBl) 

czu=coi*(ccTuu-os*nc  row) 

CZW=-CCl*nCTDW 
CZ  A = CZW*  7 
CZTC=-CO/*OC  fOT 
CZ  B 1 =-CC2 *DC  TDB 1 
C03=CCI /^CTR 


CMUR*C03* ( XR* < DCTOU-DS*DCTDW)+HR* (DCHDU-DS*OCHQW ) ) 

IF  I V ) 70»  7C»  71 

70  C04*0. 

C05*0. 

CYBB*-.0C5*I VOLF**. 67) /AO 
ENVF*0. 

C20*0. 

CMUF*0. 

CMUT*0. 

GOTO  71 

71  C04«HF*V0LF*RH0*V/ (DP*AD*ROTR ) 

ENVF  »-C01*HF*V0LF*2. *OMUHS/ ( AD*5. 73*SI GMA*2. *RGTOR ) 

CYBB*<  I-.20*CDA/( AD*.6)  )-2.86*SVT*RH0*V*V/t2.*0P*.V,3) ) 

CYVB*CYBR/V 

C05=-ST*TL*2.86*RH0*V/(0P*AD*R0TR) 

C20=-2.86*SVT*V*V/(2.*VT!P*VTIP*AD) 

CMUF*CO4*<(6LW/OMUHS)+2.*(BlC-0.O*DCTDU-3.*0MUHS*ICTAS-fcLW/2.) ) I 
CMUT*C05*< (ELW/CMUHS)*2.*«BlC+TI*DCTDU-3.*0MUHS*(CTAS-ELW/2.)n 
73  CONTINUE 

C CMU  FOR  V GT  20. 

CMU*CMUR  KMUF+CMUT 

C06=-W*2.86*ST*TL*TL/(VT!P*VnP*CTAS*5.73*SlGMA*AD*DP*AD*RGTR) 
CMWD*C06*( l.-2.*0CT0W) 

CMtoRaC03* I XR*DCTDW*HR*DCHDW) 

CMWT*CO  3* ( ( OMUHS/ ( ASA* AO) )*( VOLF*HF- (4. 3P*ST  *TL/2. ) ) ) *DCTDW 

CMW*CMWR*CNWT 

CMAD*C06*V<I1.-2.*DCTDW) 

CMAR*C03*V*(  XR*OC  TOW*HR*DCHOW) 

CMAT*C03*V*( ( OMUHS/ ( ASA* AO) ) * ( VOLF*HF- (4. 3Q*ST*TL/2. ) ) ) *OCTDW 
CMA*CMAR+CMAT 

CMTO«-C03*TL*2,86*0MUHS*ST*TL*PPF/l2.*AS4*AD)+C01*DCHDQ*HR/R0TR 
CM  IT*- ST* TL* 2. 86/ ( AO*ROTR ) 

CMIT»CMIT*RHC*V*V/(2.*0P) 

C06=W/CCTAS*OP*AD*R0TR) 

CMTC*COA*I DCTDT*XR*HR*OCHDT) 

CMB1 -C06* ( XR+OCTDR i*HR*DCHDBl ) 

C END  LONGITUDINAL  DERIVATIVES  START  LATERAL 
Y01=1.55*CMUHS**.8l 

Y0P=Y01-(10.*ELS)*(Y01*.82*0MUHS**l.R2)/12. 

Y02=-l.5*2.*0MUHS 

Y00P=Y02-U0.*ELS)*(Y02-1.5-.05*0MUHS)/12. 

Y02=.90*CMUHS**.68 

YIP=Y02-(10.*ELS)*(Y02MI. 85* OMUHS* *2. 7 ) ) / 12. 

Y2P=-.23ft*0MUHS-< 10. *5LS )♦ I . 108*0MUHS) / 12. 

Y02=-1.5-».  5*  OMUHS 

Y3P=Y02-1 10. ♦ELS)*(Y02-l.+3.38*0MUHS**1.93  )/!•). 

Y02=22.5-13.*0MUHS 

Y4P=Y02-(!0.*ELS)*(Y02-5.)/l2. 

AIC=BIC 

0CYDT=DCT0T*(  AlC-*ELi*«l.  + AS4*0Ll)CT/2.  ) ) 
0CY0V*-.0125/it.46-Y4P*CTAS*CTAS-DCTDW*AlC*< AlC*i.5*ELl ) 

DCYDT=PC YDT+DCTOT*HMF*CTAS*YQP*  <Y1P*2.*AS4*Y2P*DLDCT) 
DCYDP=.0625*CMUHS*ELl*0CT0P*(AlC*ELl*(i.*AS4*DLDCT/2.  ) ) 
OCYPP=OCYOP-I6.*CTAS*Y3P/BMF*DCTDP*BMF*(YIP*2.*AS4*Y2P*DCLDT)*CTAS 
DCYDV=DCYDV-0CTnw*AlC*BMF*CTAS*(YlP»?.*Y2P) 

CYVR=CJi*DCYOV 
CTWTR=.2S/( 1.-AS5*.25) 

CYVTR=COI*CTV«TR*SIGMAR*TRR*TRR/(SIGMA*ROTR*ROTR) 

CYVTR=-C  Y VTR 

CYB  = V*(CYVR+CYVTR ) +C  YRB 
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CYAIC=W/!PP*AD) 

CLB»C03*V*HR*DCYDV 
CLTD=COl*HR*OCYDP 
CLA1*W*HR  Z( DP*AD*ROTR ) 

ENVTR*C03*TL*CTWTR*5IGHAROTRR*TRR/(SIGMA*R0TR*R0TR) 

ENVT  *£NVF*1  • 43*SVT*TL/  !-*HF*VOLF ) 

ENVR»C03*0S*HR*DC YDV 
CN8=V*<ENVF*ENVT*ENVTR-*-ENVR> 

ENRTR*-TL«ENVTR 
ENRT  »-TL*ENVT 

ENRRR*ROTR*OS*OS*HR*OC YDP*C03 

CNSO»ENRT+ENRTR+ENRRR 

CNTD«£NRRR/DS 

CNAl»CLAl*OS 

CTTTR«l./(3.*(l.-AS5*.25n 

CNTTR»-W*TRR*TRR*TL*CTTTR/!CTAS*ROTR*ROTR*DP*AD*ROTR) 

FORMAT! 7F 10.4 ) 

CYV* ( CYVR^CYVTR ) ♦CYVB 
CLV=C03*HR*0CY0V 
CNV=ENVF +ENVT+ENVTR+ENVR 
CYSD^-CYVTR*  TL+C20 
WRITE16, Jl 

3 FORMAT ( l AX, l IHOUTPUT  DATA) 

CL$0»O. 

CYT0«O. 

WR1TE!6,4> 

4 FORMAT! 1CX, 71HL0NG! TUOINAL  DERIVATIVES  ARE  IN  ORDER  CXU,CXA,CZU,CZ 
1 A*  CMUf CMA»CMAO»CM  TD ) 

WR I TE ( 6, 5 ICXUtCXA tCZUtCZA tCMUtCMAtCMAD tCMTD 
•i  F0RMAT!R:i4.7) 

WRITE!6»6)CZttl,CMBl,CX8l 

6 FORMAT! 10X,19HCYCLIC  C0NTR0LCZBi=,ElS.8,lX,SHCMRl=,E15.8, IX.SHCXBl 
1*,E15.8) 

MRITE!6,  MCZICtCMTCtCXTC 

7 FORMAT! I0X.23HC0LLECTIVE  CONTROLCZTC*,El5.8, IX.5HCMTC*. F IS. 8, IX,5H 
1CXTC*,C1S.8» 

WRITE(6,8ICMIT 

8 FORMAT!  l.)X,25HH0RIZ.  TAIL  CONTROL  CMIT=,E15.8) 

WR ITE ( 6, 7 ) 

9 FORMAT! I JX.84HLATERAL  STABILITY  DERIVATIVES  IN  ORDER  ARC  CYB,CLR,C 
INB  tCYSDtCLSDtCNSD  tCY  TO  »CLTD  »CNTD) 

WRITE!6,5)CYR,CL&,CNR,CYS0,CLS0,CNSD,CYTD,CLT0,CNTD 
WRITE! 6, IC)CYA1C,CLAI,CNAI 

10  FORMAT! 10X,20HCYCLIC  CONTROL  CYAl« tF15. 8 , IX, 5HCL Al*f E15 .8, IX,5HCNA 
11*, 615. 01 

WRITE!6,  IDCNTTR 

11  FORMAT! 10X,28HTAIL  ROTOR  COLLECTIVE  CNTTR* ,E 15. 8 ) 

WRITE16.80) 

WR  ITE  ( 6,  •>  1DL0U , DLDW,  DLOC  T tDCTDT ,DCTDW, OCTDU , DCT DBl  , DCTDP 
WR l TE ( 6, R 1 1 

WRITE(6,S)0CHDT,DCHDW,DCHDU,DCHDBl 

WRITli(6,82) 

BIC  = BIC*  j7.  i 

WRITE!6,c»)fUC,CTAS,CTfELOfEU,FLS,CQOS,TlC 
RO  FORMAT!  IX, 68HAUXILIARY  DERIVATIVES  DLOU  , DLOW  , OL  OCT  , 0CTDT  , DCT  f'W , DC T 
10U, OCTDB l » DC !DP  ) 

PI  FORMAT! IX, ?5H  DCHO T» OCHDW , DCHDUtDCHDBl ) 

82  FORMAT!  IX, 3 B 1C ,CT AS,CT ,ELO,ELl , ELS, COOS , T 1C ) 

WR ITE ! 6, 76 ) 

WRITE!6,c')CXW,CZVi,CMW,CMW0,CYV,CLV,f.NV 
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96  FORMAT! 10X,69HTHE  VELOCITY  DERIVATES  ARE  IN  ORDER  CXW,CZH*CMW,CMWD 
1*CYV»CL V,CNV ) 

93  FORMAT! 20X» IflH  HELICOPTER  SYSTEM} 

99  FORMAT! 20X, 16H  AIRCRAFT  SYSTEM) 

95  FORMAT! 2 5X»5HH0VER ) 

50  FORMAT! IHC) 

WR ITE! 6*  5)CMUR  »CMAR ,CMWR  »CZH»CZA 
HR  IT E ( 6«  50) 

IF1VU9,  15,10 
19  CONTINUE 
CLTTR.O. 

CYTTR*-R*CNTTR/TL 
C STATIC  STABILITY 

C STEAOY  SIDESLIPS 

DTOB«-tOP*AD/W)*!CYB*<CLTTR*CNAi-CNTTR*CLAl)-CLR*!CYTTR*CNAl-CNTTR 

1*CYA1C)*CNB*!CYTTR*CLAI-CLTTR*CYAIC)) 

DTROT*--!  CLB*CNA1-CNB*CLA1 ) /DTOB 
OAIDT«-!CLB*CNTTR+CNB*CLTTR)/OTDB 
DN2=CLTT«*CNA1-CNTTR*CLA1 
DTOB-OTDB/DN2 

HR ITE !6« 18)DTDB«DTRDT«0AIDT 
C STEAOY  TURN 

PHI*0« 

DPHIM0./57.3 
DO  15  K*l,5 
PHI -PHI+DPHI 

SN*32«2*SIN(PHl ) / (COS! PHI )*V) 

THT ■- 30. /57. 3 
DTHT-15./57.3 
00  15  KK*1, 3 
THT«THT*DTHT 

TTR*-SN*!THT*(CLA1*CLT&-CNAI*CNTD)-C0S(PHI )*<CLA1*CLSD-CNAI*CNSD> ) 
TTR*TTR/0N2 

Ai*-SN*( THT* !-CLTTR*CNTD*CNTTR*CLTD) -COS (PHI )*(  CLTTR*CNSD*CNTTR*C.L 
ISO ) »/0N2 
PF»1./C0S(PHI) 

Bl«SN*SIN(PH! )-»( ! 1 1.  /COSIPHI ) )-l.  ) *CMA*W/  IOP*AD ) ) / i U A*CMTD-CMA*CZ 
1T0 ) 

THP=THT*57«  3 

16  FORMAT! 1X« 12H  STEAOY  TURN.12H  BANK  ANGLF* ,F5 .0, l 3H  CLIMB  ANGLE*?  F5 
1.0.13H  L 'AO  FACTOR*, F6. 2) 

PHP=PHl*57.3 

MR ITE 16? 16)PHP»THP»PF 

TTRP*57«.  3*TTR 

AlP=Al*57.3 

BIP=B1*57.3 

WRITE(6,  17)TTRP,AIP,B1P 

17  FORMAT! IX»38H  CONTROL  ANGLES  TAIL  ROTOR  COLLECT IVE*?F6. 2, 16H  LATER 
1AL  CYCLIC*, F6. 2, 1AH  PITCH  CYCLIC*, F6. 2) 

18  FORMAT! IX, 33H  STEAOY  SIDESLIP  RATIOS  PHI/BETA*, F9. 5, 27H  TAIL  ROTOR 
1 C0LLECT!VE/PHI=,F9.5,2OH  LATERAL  CYCLIC/PHI  * ,F9. 5 ) 

15  CONTINU= 

G0T077 

END 

SUBROUTINE  TANNRH ( TI PMH,  S I GMA ,CTOY.H,CQOSH,CTSTH,  IREAD) 

DIMENSION  SIM  7 ) ,CTH ( 9, 5 , 3 > ,CQH<9,5 ,3) ,CQX ( 3 ) ,CQXX ( 3 ) 

EQUIVALENCE! IN, INPUT) 

150  FORMAT! 13F6. A) 

IN  = 5 
MADLY=6 
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IF!IREAD)153,153,154 

153  READ!  INPUT, 150) ! ! <CTH!K, J, I ) ,K*l,9) , J=l ,5 ) , I =1 » 3 ) 

READ( INPUT, 150) < ( (CQH! K*  J* I 1 ,K*1 ,91 , J*l,5) ,1*1,3) 

SIG( i )*.05 

00  151  1*2,6 
151  SIG!I)*SIGU-1K.C25 

149  FORMAT! 1H  , 38HTANNER  DATA  FOR  HOVER  HAS  BEEN  READ  IN) 

WR I TE! MADLY, 149) 

GO  TO  50 

154  00  155  J*l,4 

IF ( S IG<  J )-S IGMA )155,156,156 

155  CONTINUF 

156  CALL  ZINDEX! J,JA, JB, JC) 

00  157  K *1, 8 

IF(CTH(K, l, 1)-C TOSH >157,158,158 

157  CONTINUE 

158  CALL  ZINOEX(K,KA»KB,KC) 

00  160  1*1,3 

00  159  J=JC,JA 

CALL  CURVEB(CTH!KA,J,I ) ,CTH(KB,J,I ) ,CTH!KC , J , I ) ,COH(KA, J,  I ) , 
1CQH(KB«J,I),CQH(KC*J,!)*CY0SH,CQX(J)) 

159  CONTINUE 

CALL  CUR VEB (SIG(JC),SIG(JR) ,SIGI JA) ,COX( JC) ,CQX< JB) ,CQX( J A), SIGMA, 
lCQXX(D) 

160  CONTINUE 

CALL  CURVES! G, 5, G, 6, 0«  7,CQXX (1) ,C0XX(2),CQXX(3) ,TlPMH,CQOSH) 

CALL  CURVEBtSIG! JC),SIG! Jft),SIG! JA),CTH!9,JC,l) ,CTH!9,JB,l), 
1CTH!9,JA,D, SIGMA, CTHX) 

CALL  CURVES (SIG(JC),SIG(JB),SIG(JA) ,CTH!8, JC,2) ,CTH!«,JB,2), 
ICTH!8,JA,2),SIGMA,CTHY) 

CALL  CURVEB(SIG!JC),SIG(JB),SIG!JA) ,CTH(7,JC,3) ,CTH(7 , JB, 3 ) , 
ICTH!7,JA,3),SIGMA,CTHZ) 

CALL  CURVEB!C.5,0.6.0.7,CTHX,CTHY,CTHZ,TIPMH,CTSTH) 

CQOSH*l.l*CCKSH 
50  RETURN 
END 

SUBROUTINE  CURVEDIX1,X2,X3,Y1,Y2,Y3,X,Y) 

IF!Xl-X2)iH,llO,lll 

110  XWU.OOOl 
HRITE(6,120)XI,X2,X3,Y1,Y2,Y3 

111  IF(Xl-X3)U3,ll2,ll3 

112  Xl=X1^.0CCl 
WRITE!6,l20)Xl,X2,X3,Yl,Y2,Y3 

113  IFK2-X3!  115,114,115 

114  X2**2*.CCCl 

WRITE (6, 120)Xl,X2,X3,Yi,Y2,Y3 

115  CONTINUF 

120  FORMAT! 6F 15.8) 

Bl*Yl/( !Xl-X2)*!Xl-X3) ) 

B2*Y2/! (X2-X1 )*!X2-X3) ) 

B3=Y3/!(X3-Xl)*!X3-X2)) 

Al=!Bl*X2*X3>*!B2*Xl*X3)MB3*Xl*X2> 

A2*-B1*( X2+X3J-R2*! X 1 + X3 ) -B3* ! Xl»X2> 

A3=B1*R’«P3 

Y*A1»(A2*X  ) MA3*(  X**2)  ) 

RETURN 

END 

SUBPOUTI  iE  r.U9V-C!Xl,X2,X3,Yi,Y2,Y3,X,Y) 

Z0=( (Y2-Y1 ) /! XI-X2) )-( !Y3-Yl )/! XI -X3) ) 

IF!  ABSIZD)-...  1C00001  11,2,2 
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1 Z0«. 00000001 

2 ZA»(  (((Vl*Xl)-<Y3*X3n/(Xl-X3n-(((Yl*Xl)-(Y2*X2)  »/<Xl-X2>  n/zo 
ZC*( (YlOXl)-!  Y1*ZA)~(Y3*X3)'MY3*ZA) )/(Xl-X3) 

Y«zc*tnxi-zA)*(Yi-zcn/(x-zAn 

RETURN 

END 

SUBHOUTINE  ZINDEX(!T,ITA,ITB,ITC) 

ITA-IT 
ITB-IT-l 
ITC*IT-2 
IFUTCJltl*2 
l ITA*3 
1TB«2 
ITC*1 

2 CONTINUE 
RETURN 
cNC 

SUBROUTINE  TANNSS < TI PM,0MU,CL ,CD,CQ,CDM7, l RE AD, AC ) 

DIMENSION  CDS(4,6,25),CQS(4,6,25)  ,EAC(4,6,25) , Jl  < A ) ,PU(  A ) 

DIMENSION  CL  A (4,6, 5)  ,CDM(4,6,5)  tK|  (5)  ,C07(5»  ,C07  (5)  ,EA7<  5 ) ,PCC!  5 ) 
BOE  ROTOR  PERF.  CHARTS  STRAIGHT  LINF  INTERPOLATION  VALUES  OF  MT/MU 
MT«.6fl/.2,.5,.8,l.,l.2,t.5  MT-.95/. 31 ,.47, .66 

MT«. 77/. 33, .52,. 79, I. 06  MT». 875/. 35, .52, .75, 1.05  VARIOUS  CL 
INTERPOLATE  BETWEEN  MT  IS  DM,  OMU  IS  PU,  CL  IS  PCC 
NOTE  IF  BLADE  IS  STALLED  ALPH  14  OR  Al  SIGMA  IS  INCREASED 
NOTE  IF  OMU  IS  HIGH  ROTOR  THRUST  IS  LOW  PROP  THRUST  IS  ADDED 
IN  READ  (Nl,J,K)  Nl  IS  MT,  J IS  OMU,  K S VALUES  OF  ALPH  FOR  EACH 
OF  5 GROUPS  OF  CL  CONTROLLED  BY  JJ 
I F < I READ- 1 1 86,70,  70 
86  NI«1 
N2«6 

1 D02J-I.N2 

8 7 FORMAT!  I'-' !5> 

!A»1 

!B«5 

002J J»l, 5 

READ(5,104HCDS(NI,J,K),K>IA,IB) 

REAn<5,104)(C0S(Nl,J,K),K*IA,IB> 

REAO(5,104)(EAC(NI,J,K),K«IA,1B) 

READ(5»104)CIA(N1,J,JJ),CDM(N1»J»JJ) 

I A* I A+5 
IB»IB>5 

2 CONTINUE 
N 1 *M 1 ♦ 1 

IF (N 1-31 3, 3,4 

3 N2*4 
GOTOI 

4 IFIN1-4) 5,5,6 

5 N2* 3 
GOTOI 

6 WR ITE(6»91 ) 

G0T090 

91  FORMAT! 1H  ,11HB0E  DATA  IN) 

70  N l = 1 

I F C T I PM- .77)7,7,8 

7 11  = 1 

PM=(TIPM-.67S)/.095 
GOTO! 1 

8 IFHTIPM/.87S  )-l . )9,9,10 


PM=(TIPM-.77)/.l05 

G0T025 

10  11*3 

PM*(TlPM-.875)/.C7‘.> 

G0T035 

11  IF1PMU- .5)12,12, 13 

12  J UNI  1*1 

PU1N1  ) = ( 0MU- . 2)  / . 3 
GOT020 

13  IF10MU-. 8)14, 14,15 

14  JKNl  )=2 

PUIN1 1 *<  OMU-. 5 1 /. 3 
GOT02O 

15  IF10MU-1. 016,16,  17 

16  JKNl  )=3 
PU(Ni)*OMU-.8)/.2 
G0T020 

17  |F(0MU-1. 2)18,18,19 
Ifl  Jimi=4 

PU(N1)*( . MU-l.C  ) / • 2 
GOT020 

19  J 1 1 N 1 ) = 5 

PUINi  )*< : HU-1.21/. 3 

20  X2«2 

N 1 « N 1 ♦ l 

25  IF(0MU-.52)26,26,27 

26  JICN1 >*l 
PUni)*(CMU-.33)/.l9 
G0T030 

27  1F(0MU“. 79)28,28,29 

28  J 1 ( Ml  1=2 

PU  ( N 1 ) = ( ■''HU-  . 52  1 / . 2 7 
G0TC30 

29  J11N1)»3 
PU(Nl)*<v>U-.79)/.37 

30  IF(Nl-2132tlCO,lOC 
32  X2*3 

N1*NU1 

35  IFOHU-. 52)36, 36,37 

36  Jl(Nll=l 

PU  ( N 1 1*1 OMU- . 35 1 / . I 7 
GOT  340 

37  IFICMU-. 7513ft, 38, 39 

38  J1 (N1 )*2 

PUINI  1*1  ')HU“.52)/.23 
G0TC40 

39  J1(M 11=3 

PUCNl >=<  1HU-.751/.3 

40  IF1M1-2142, 100, IOC 
42  X2  = 4 

N l =*!  l-*l 

42  IFIOMU- .47)43,43,44 

43  Jl ( N1 1*1 

PU(N1)=( JMU- . 32 )/• 15 
GO  TO  IOC 

44  Jl INI ) *2 

PU  (N 1 ) = ( OMtJ-  . 47  ) / . 19 
iro  J= J 1 c l ) 
i = il 

0054KK=1 . 4 
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D045Kl»l,5 

N*Kl 

IF (CLA(I,J,KU-CL  >45,45,46 

45  CONTINUE 

46  IF(N-II47,47,48 

47  N»2 

48  DEN=CLA<  I,J,N)-CLAU,J,N-1) 

PCC(KK)»(CL-CLAI I , J, (N-l ) ) l/OEN 
K 1 (KK  )«i*5*N-5 
IF(KK-l)49,49,5) 

49  J»JI(  !)♦! 

GOT054 

5C  1F(KK-2>51,51,52 

51  l«IU 
J-J1(2) 

GOT354 

52  IF<KK-3)53, 53,54 

53  J*Jl ( 2)4 1 
44  CONTINUE 

U2*im 

0055  KK«t,5 

Kil»Kl<l)*KK-6 

K 12*K 11*5 

K 1 3*K  1(  2 J4KK-6 

K14=K13*5 

Kl5=Kl(3)4KK-6 

K16*K 15*5 

K17*Kl(4MKK-6 

K18-K1745 

JX»Jifl) 

JY»J 1 ( 2) 

JW»JX*1 
JV*  JY ♦ 1 

CDl*COSIIl,JX,Kll)*(l.-FCCIU)4PCCm*COSUl,JX,K12) 
CD2*C0SIIl,JW,Kl3)*U.-PCC<2>)4PCCt2)*CD$< U,JM,K14) 
C03*C0S(  112,JY,K15)*(l.-PCC(3))4PCCI3)*C0S(U2,  JY,K16) 
C04»C0S(  I12,JV,K17)*(1.-PCC(4)WPCC(4»4C0S(  112,  JV,K18> 
C05«C0l4PUt  1 )M  CD2-C01 1 
CD6=C03*PU( 2 )♦ (C04-CD3 ) 

C07(KK)»C054PM*(C06-f.D5) 

CQ1=C0S(  Il,JX,Kin*(l.-PCC(l)  )4PCCU)*C0S(U,JX,K12) 
CQ2«C0Sm,JW,Ki3)*Il.-PCCI2))4PCC(2)*C0Sm,JW,Kl4) 
CQ3«C0SI  U2,JY,K15)*II.-PCCI3))4PCCI3)*CQSIU2,  JY,K16) 
CQ4*C0S( 1 12, JY,K1 7 )*l l.-PCC  f 4) ) ♦PCCI41 4CQSI 1 12, JV,K18 ) 
CQ5»CQl4PUIl)*(CQ2-CQl) 

CQ6=C034PU(2)*IC04-CQ3) 

CQ7(KK»»CQ5*PM*IC06-C05) 

EACl*EACUtvJXvKll)*f  l.-PCC(l)  )4PCC«l)*«:ACUl,JX,Kl2) 
EAC2=EACf Il,JW,Xl3)*f 1.-PCCI2) )4PCCI2)*EACf I1,JW,K14) 
EAC3*EACII12,JY,K15)*(1.-PCC(3»)*PCC13)*FACU12,JY,K16) 
EAC4«EAC< I12,JV,K17)*I1.-PCCI4))4PCCI4)*FACI 112,JV,«18) 
£A5=EACl»PUm*I£AC?-EACl> 

EA6=EAC3*PU( 2)*(EAC4-EAC3) 

55  EA7(KK)=; A5*PM*< EA6-EA5) 

N 1 1* ( K If 1 1 *4  ) /5 
N12=N).  1-  l 
N13=(KH2>*4)/5 
N 14=N  13-1 
N15=(KH  3H4I/5 
N16=N15-  1 
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N17=(KU4)*4)/5 

N18*N17-1 

coi*CDM(n,jx,Ni2>*(i.-Pcc(i»H-PCC(i>*coM{ii,jx,Nm 
CQ2=CDM( I1,JW»N14)*( l.-PCC(2 ) l+PCC (2 ) *CDM( 1 1 , JW ,N13) 
CD3*COM(I12,  JY,Ni6)*(l.-PCC(3)  > +PCC  (3> *CDM( 1 12,  JY,N15) 
C04*C0M( 112, JV,N18»*( 1.-PCC14H *800(41 *COM( 112, JV,N17) 
C05«C01*PU( 1 >♦(002-00! > 

CD6*CD3*PU(2)*(C04-CD3) 

C0M7*C05*PM*(C06-CD5 ) 

D057KK*1, 5 
N-KK 

DEN*CD-CD7(KK ) 

IF(C0-C07(KK) >57, 157,157 

57  CONTINUE 
G0TC60 

157  1F(N-1)5H,58,59 

58  N«2 

59  P* (CO-CO 7 (N- il)/(C07(N)-CD7(N-l)> 
AC«EA7n-l)*P*(EA?(N)-EA7(N-l>> 

C0»C07(N- 1 )*P* (CQ7 (N )-CQ7(N-l ) ) 

G0TC89 

60  AC-10000. 

C AC-10000.  FOR  ROTOR  CONDITION  Of  INADEQUATE  THRUST 

89  CONTINUE 

88  FORMATI 3X7E15.8) 

K4  FORMATI 10F13.51 

90  RETURN 
END 

* 1 1 ST( STOP  ) 
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.OiC 

.0126  .0174 

.00055.002  .0034 

.0051 

.0070 

.0090 

.0116 

.015 

.0190 

.00055.0921 

.0037  .0055 

.0076  .0104  .0132 

.317 

.0204 

.00055 

>.0922 

.3  >4 

.006 

.CC82 

Oil 

.0146  .0189 

.0215  .00055.00162 

.00275 

.004 

.0053 

.0071 

.0092 

.012 

.0134 

.0C055 

.0018  .0031 

.0045  .0064  .0085 

.Oil 

.0145 

.0155 

.00055 

.0020 

. OC  34 

.0051 

.0071 

.0084  .0125 

.0144  .0168  .0055 

.0021 

.0037 

.0055 

.0078 

.0135 

.:i4 

.016 

.018 

.00055.0022 

.004  .006  .0085 

.0115 

.0154 

.018 

.01871 

.0055 

•CC165 

.0028 

.0042 

.0049  .00575 

.00665.00775.0104 

.00055 

>•0016 

.0032 

.0049 

.0058 

.00675 

.0079 

.0  .'92 

.0131  .00055 

1.00195. 0035  .00525 

.0064 

.0076 

.0088 

.0105 

.0140 

-C :55 

.0021 

• Qw  ? 8 

.0057  .0068 

.0084  .0098  .0115 

.9148 

.00055 

.0022 

.0040 

.0060 

.00725 

.0088 

.0106 

.0127  .0154 

.00(25 

• 

0002 

-.001 

-.0025  -.0045 

l 

.00084 

• 

091  9 

.0.149 

.00183  .00233 

10. 

0 

• 

-10. 

-20.  -30. 

.01 

- 

8<~. 

.00575 

• 

001 

-.003 

-.00825  -.0145 

2 

.0001 

• 

0014 

.00203  • 

.00307  .00441 

10. 

0 

• 

1 0 . 

-20.  -30. 

.01 

.01175 

-.0003 

10. 

.05 

.011*25 

-.00084 

10. 

.07 

.021 

-.00109 

10. 

.09 

.00131 

.00117 

10. 

.01 

.00311 

.0005 

1C. 

.02 

.00507 

-.00042 

10. 

.03 

.00703 
-«0C 134 
1C. 

.04 

. 0^916 
-.00334 
10. 

.05 

.0032 

.00132 

12. 

.01 

.0041 

.00079 

12. 

.02 

.00525 
.CGI  32 
12. 

.03 

.00701 

-.00106 

12. 

.04 

.01189 

-.00201 

12. 

.05 

.00597 

.00051 

12. 

.01 

.00698 
- 127 
12. 


-60. 

.002 

.00119 

0. 

-60. 

.003 

.00139 

0. 

-50. 

.00575 

.00183 

0. 

-33r 

.00082 

.00142 

0. 

-70. 

.00147 

.00109 

0. 

-55. 

0.00213 

.00075 

0. 

-40. 

.00327 

.000307 

0. 

-31. 

•C04C9 

.00008 

0. 

-25. 

.00205 

.00148 

fi. 

-17.5 

.00262 

.00119 

8. 

-15.8 

.00361 

.00059 

8. 

-10.5 

.00513 

-.0004 

8. 

-11. 

.00656 

-.00125 

8. 

-9. 

.00343 

.00127 

8. 

-14. 

.00^13 

.OG-^ 

8. 


-.006 

.00277 

-10. 

-.008 

.00376 

-10. 

-.00925 

.0047 

-10. 

.00033 

.00176 

-10. 

-.00048 

.00209 

-1C*. 

-.00098 

.00242 

-10. 

-.00131 

.00268 

-10. 

-.00131 

.00268 

-10. 

.00184 

.00158 

4. 

.00205 

.00142 

4. 

.00262 

.00106 

4. 

.00389 

.00023 

4. 

.0049? 

-.00053 

4. 

.00267 

.00140 

4. 

.00305 

.00140 

4. 


-.01475 

.00455 

-20. 

-.01975 

.00614 

-20. 

-.024 

.00817 

-20. 

0. 

.00209 

-20. 

-.00245 

.00334 

-20. 

-.004780 

.00451 

-20. 

-.00638 

.00527 

-20. 

-.00752 

.00602 

-20. 

.00184 

.00158 

0. 

.00205 

.00135 

0. 

.00258 

.00106 

0. 

.00348 

.00026 

0. 

.00438 

-.0003 

0. 

.00267 

.00178 

G. 

.00305 

.00127 

0. 


-.02425 

.00658 

-30. 

-.03325 

.00921 

-30. 

-.0345 

.01064 

-27. 

-.00131 

.00318 

-30. 

-.00523 

.00502 

-30. 

-.00948 
. C0736 
-30. 

-.01226 

.00903 

-30. 

-.0104 

.00752 

-25. 

.00226 

.00158 

-4. 

.00238 

.00132 

-4. 

.00246 

.00116 

-4. 

.00328 

.00053 

-A. 

.00402 

.00023 

— A, 

.00337 

.00133 

-4. 

.00330 

.00051 

-4. 


.03 

.00921 

-9.5 

.00603 

.00419 

.00413 

-.oom 

-.00025 

.00006 

-.00025 

12. 

8. 

4. 

0. 

.05 

.0127 

-6. 

.00889 

.00698 

.00698 

-.00381 

-.00241 

-.00152 

-.00127 

12. 

8. 

4. 

0. 

.07 

.01778 

-3.3 

.0141 

.01111 

.0111 

-.0  0 711 

-.00521 

-.0033 

-.00127 

12. 

8. 

4. 

0. 

.09 

.00712 

-.5 

.00398 

•0v35l 

.00527 

0. 

.00495 

.00165 

.00088 

ft. 

4. 

0. 

-4. 

.01 

.C 1027 

-12. 
.0048 l 

.00370 

.00416 

0. 

.00132 

.00165 

0. 

8. 

4. 

0. 

-4. 

.0  4 

.01341 

-7. 

.00666 

.00388 

.00462 

-.00319 

0. 

.00132 

.00066 

12. 

8. 

4. 

0. 

.05 

.01554 

-5.5 

.00897 

.0C555 

.00065 

-.00429 

-.00110 

.C0O22 

-.0011 

12. 

ft. 

4. 

0. 

.07 

.01536 

-3.7 

.01221 

.0086 

.00999 

-.0066 

-.00341 

-.00209 

-.00165 

1C. 

8. 

4. 

0. 

.09 

• C .77 

-2. 

.00598 

.0057 

.00656 

.0002 

.00064 

.0015 

.00086 

4. 

2. 

0. 

-2. 

.01 

.00627 

-8. 

.00513 

.00542 

.00627 

.0002 

.0015 

.00193 

.00064 

4. 

2. 

0. 

-2. 

.03 

.01083 

-5.3 

.00593 

.00513 

.00556 

-.00129 

.005 .3 

.00193 

.00129 

H. 

4. 

2. 

0. 

.05 

.Cl  111 

-4  3 

.01198 

.0057 

.00684 

-.  >'193 

.o:)64 

.0017? 

0. 

ft. 

4. 

2. 

0. 

.07 

.01  III 

-3  7 
.G/61-* 

.0.656 

.00855 

- .OC.’  79 

.0064 

.00107 

-.0017? 

8. 

4 

2. 

0. 

.09 

.00173 

-.‘.6 

.C002» 

-.00035 

.00139 

.00091 

0014 

. OC l 37 

.00182 

10. 

J. 

-lo. 

-20. 

.01 

.00554 

-80. 
.001 39 

- . 0C2  7 7 

-.00693 

-.00023 

.OOK  3 

.00239 

.00376 

10. 

0. 

-l :. 

-20. 

A 1 

14 


.00540 

0. 

-4. 

.0087 

-.0002 

-4. 

.0139 

.0002 

-4. 

.00025 

-.00132 

-P. 

.00601 

-.0022 

-B. 

.00555 

-.00055 

-2. 

.00786 

-.00165 

1. 

.0125ft 

-.00143 

.00883 

-.00021 

-4. 

.00926 

-.00064 

-4. 

.0081? 

-.000143 

-2. 

.001 

-.0017? 

-2. 

.01126 

-.00215 

l. 

.00346 

.00274 

-3">. 

-.01247 

.00593 

-30. 


.03 

•C1C39 

-.00068 

10. 

.05 

.01596 

-.00174 

10. 

.07 

.01951 

-.00228 

5. 

.08 

.00177 
.00 117 
10. 

.01 

.00319 

.0004 

10. 

.02 

.00478 

.00070 

10. 

.03 

.0069 

-.0014 

10. 

.04 

.0089 

-.00182 

5. 

.05 

.00211 

.00159 

8. 

.01 

• Of  235 
.00122 
8. 

.02 

.00541 

•0CC06 

12. 

.03 

.00689 

-.00083 

12. 

.04 

.00864 

-.00184 

12. 

.05 

.0:533 

.00105 

8. 

.01 

.00568 

.00084 

8. 


-58. 

. 00:98 
.00068 
0. 

-45. 

.00589 

.00023 

0. 

-28. 

.00728 

-.0004 

0. 

-20. 

.00133 

.0014 

0. 

-62. 

.00169 

.00112 

0. 

-48. 

.00195 

.00084 

0. 

-39. 

.00301 

.00037 

0. 

-31. 

.00407 

-.00009 

0. 

-23. 
.001 79 
.00159 
0. 

-16. 

.00223 

.00135 

0. 

-12. 

.00287 

.00092 

5. 

-12. 

.00398 

.00024 

5. 

-11. 
.00549 
-. 00*61 
5. 

-9.5 

.00361 

.00136 

4. 

-13. 

.00473 

.90126 

4. 


-.00416 

.00235 

10. 

-•0C45 

.00353 

-10. 

-.0045 

.00388 

-10. 

.00009 

.00164 

-10. 

0. 

.00187 

-10. 

-.00071 

.00224 

-10. 

-.00089 

.00252 

-10. 

-.00106 

.0-262 

-10. 

.00239 

.00147 

-5. 

.10219 

.00153 

-5. 

.00271 
.CO  101 
0. 

.00342 

.00043 

0. 

.00438 

-.00024 

0. 

.00344 

.0021 

0. 

.00344 
.00168 
0 . 


-.01247 

.00559 

20. 

-.01559 

.00752 

-20. 

-.01745 

.00906 

-15. 

.00018 

.00206 

-20. 

-.00212 

.00318 

-20. 

-.00372 

.0043 

-20. 

-.00531 

.00514 

-20. 

-.00655 

.00608 

-20. 

.00386 

.00098 

-10. 

.00318 

.00129 

-10. 

.00239 

.00138 

-5. 

.0033 

.00098 

-5. 

.00386 

.00061 

-5. 

.00568 

.00084 

-5. 

.00559 

.00084 

-5. 


-.0201 

.00866 

30. 

-.02772 

.C1254 

-30. 

-.02979 

.01444 

-?0. 

-.00106 

.00309 

-30. 

-.0046 

.00486 

-30. 

-.00743 

.00655 

-30. 

-.01062 

.00832 

-30. 

-.01292 

."1966 

-30. 

.00478 

.00098 

-15. 

.00386 

.00092 

-15. 

."0316 

."0138 

-10. 

.00318 

.00141 

-10. 

.00366 

.00181 

-10. 

.01178 

-.00158 

-10. 

• "1161 
-.00147 
-10. 


A-  15 


.03 

.0056R 

-.GC1U 

fi. 

.05 

•01 144 

-.00252 

1C. 

.07 

.C142R 
— . C C 3 i 
10. 

.09 

.0C 19R 
.00111 
1C. 

.01 

.0<  592 
. V - 2 

10. 

• 03 

.01027 

-.00139 

1C. 

.05 

.01303 
- . OC 195 
10. 

.06 

.01119 

-.0°)06 

10. 

.07 
.0023 
.0  181 
10. 

.01 

.OC  3 7? 
.Of  109 
10. 

.02 

.0  549 
. CCC2 
K . 

.0  3 

.00766 

-.00076 

10. 

.04 

0.<  C991 
-.00199 
10. 

.06 
.0048 
.C  413 
12. 

.01 

.0  568 
.Ci 335 

12. 


-8. 

.00378 

.00147 

4. 

-5.7 

.00602 

.00021 

5. 

-4. 

.00851 

-.00168 

5. 

-2. 

.00103 

.00139 

0. 

-RC. 

.001)8 

.001)9 

0. 

-70. 

.0034 

.00106 

0. 

-55. 

.00442 

.00080 

0. 

-47. 

.00612 

.00056 

0. 

-25. 

.00159 

.00199 

0. 

-57. 

.00186 

.00181 

0. 

-48. 

0.00266 

.00143 

0. 

-39. 

.00275 

.0C095 

0. 

-31. 

0.00478 

.00028 

0. 

-23. 

.oom 

.00409 

6. 

-19. 

. 00  3 H 2 
.00337 

6. 


.00499 

.00042 

0. 

•C00662 

-.00063 

0. 

.00791 

-.0C168 

3. 

.00039 

.00167 

-10. 

-.00047 
.0.'  278 

-10. 

-.00348 
• C0334 

-10. 

-.00395 

.00400 

-10. 

-.00395 
. CC-403 
-1.,. 

.00124 
• Ov.238 

-10. 

. 10053 
.0C23« 
-1C. 

0. 

.00266 

-10. 

-.00004 

.00285 

-1  «. 

- . OCC  35 
, Cu  3 14 

-1C. 

.00343 

.00370 

0. 

.00343 
. OC  348 

G. 


.00654 

.0002 

-3. 

.00877 

.00021 

-3. 

.00955 

-.00126 

0. 

-.00119 

.00222 

-20. 

-.00632 

.00417 

*20. 

-.01082 

.00612 

-20. 

-.01224 

.00723 

-20. 

-.01422 

.00834 

-20. 

.00071 

.00371 

-20. 

-.00142 

.00352 

-20. 

-.00319 

.00475 

-20. 

-.00478 

.00570 

-20. 

-.0062 

.00665 

-20. 

.00519 

.00335 

-8. 

.00392 

.00361 

-8. 


.00946 

.00084 

-6. 

."1238 

-.00179 

-6. 

.01273 

.''0126 

-3. 

-.00356 

.00334 

-30. 

-.01146 

.00612 

-30. 

-.019)5 

.00973 

-30. 

-.02291 

.0118? 

-30. 

-.02528 

.01557 

-30. 

-.00057 

.00541 

-30. 

-.00407 

.0054? 

-30. 

-.00726 

.00741 

-30. 

-.01027 

.00931 

-30. 

-.0094 

.0095 

-25. 

.01 

.00279 

-?0. 

.00598 

.00)70 

-?0. 


4-16 


.02 

-17.8 

.00686 

.00451 

.00249 

.00322 

12. 

6. 

.01 

-16. 

.0098 

.00666 

.00034 

.00163 

12. 

6. 

.04 

-11. 

.01372 

.00931 

-.00241 

-.00043 

12. 

6. 

.09 

-9. 

.00674 

.00483 

.00283 

.00224 

10. 

6« 

.01 

-12.5 

.00994 

.00483 

.00C92 

.00276 

12. 

6. 

.03 

-8.5 

.01193 

.00625 

- . vjOO  3 7 

.002 .1 

12. 

6. 

.or> 

-5.7 

.01926 

.00923 

- .00238 

.0004 

12. 

6. 

.07 

-1. 

.0186 

.01235 

-.00447 

-.00209 

12. 

6. 

.09 

-80. 

.00307 

.00192 

.00223 

.00244 

10. 

0. 

.01 

-70. 

.00685 

.00307 

.0-091 

.002  j3 

10. 

0. 

.03 

-45. 

.00112 

-.0)197 

-.00041 

.00173 

10. 

0. 

.05 

-28. 

.01664 

.0064 

-.00162 

.00142 

1C. 

0. 

.07 

-19. 

> .01971 

.00794 

4 -.CC223 

.00158 

10. 

0. 

.08 

-68. 

- .GC335 

-.00251 

.0. 262 

.0^ 2 75 

10. 

0. 

.01 

-46. 

.0067 

.00349 

[ . 00C9S 

.00229 

f 10. 

0. 

* 


.00392 

.00392 

.00510 

.00344 

.00387 

.00434 

0. 

“fl* 

-16. 

.00534 

.00441 

. 00466 

.00224 

.00335 

.00473 

0. 

-6. 

-12. 

.0076 

.00637 

.00539 

.00077 

.0031 

.00464 

0. 

-6. 

-10. 

.00412 

.00625 

.01022 

.00305 

.00231 

.00015 

0. 

-6. 

-12. 

.00454 

.00674 

.01235 

.00276 

.00209 

.00075 

0. 

-6. 

-12. 

.00554 

.00596 

.00817 

.00224 

.00164 

.00224 

3. 

0. 

-4. 

.00795 

.0088 

.01122 

.00224 

.00596 

.00283 

3. 

0. 

-3. 

.0118 

.01257 

.01463 

-.00112 

.00127 

.00343 

3. 

0. 

-3. 

.00004 

-.00109 

-.00307 

.00325 

.00294 

.00355 

-10. 

-20. 

-30. 

-.0016 

-.0064 

-.01152 

.30345 

.00487 

.0065 

-10. 

-20. 

-30. 

-.00307 

-.0112 

-.01984 

.00416 

.007 

.01035 

-10. 

-20. 

-30. 

-.00397 

-.01472 

-.02637 

.00518 

.00954 

.01563 

-10. 

-20. 

-80. 

-.00384 

-.0096 

-.016 

.00599 

.00832 

.01157 

-10. 

-20. 

-30. 

-.00181 

-.00098 

.00056 

.00294 

.00334 

.00438 

-10. 

-20. 

-30. 

.0007 

-.00377 

- .008^9 

.00379 

.00549 

.00798 

-10. 

-20. 

A- 17 

-30. 

.C"» 

-36. 

.00893 

.00413 

-.00042 

-.0053 

-.01116 

.0004 

.00196 

.00405 

.00719 

.01033 

10. 

0. 

-10. 

-20. 

-30. 

.04 

-28. 

.01144 

.00558 

-.00042 

-.00684 

-.01395 

-.C0078 

-.00144 

.00432 

.00778 

.01243 

10. 

0. 

-10. 

-20. 

-30. 

.C6 

-20. 

.01381 

.00698 

-.00014 

-.00405 

-.00767 

•0C 196 

.00098 

.00445 

.00654 

.00883 

10. 

0. 

-10. 

-20. 

-25. 

.06 

-52. 

.0053 

.00394 

.00316 

.00394 

.00438 

.0042 

.0051 

.0066 

.00495 

.0044 

20. 

10. 

0. 

-10. 

-20. 

.01 

-53. 

.00393 

.00349 

.00238 

.00127 

-.00055 

.00348 

.00379 

.00512 

.00659 

.00997 

0. 

-10. 

-20. 

-30. 

-4C. 

.02 

-46. 

• Of  465 

.00316 

.00122 

-.00105 

-.00377 

•0C275 

.00384 

.00604 

.00878 

.01363 

0. 

-1C. 

-20. 

-30. 

-40. 

.03 

-38. 

.00543 

.00366 

.00022 

-.00111 

-.00249 

.0'  229 

.00412 

.00714 

.00906 

.01153 

0. 

-10. 

-20. 

-25. 

-30. 

.04 

-29. 

.01.  5 

• 0 6 3 8 

.Cl332 

.00166 

-.00014 

-.00037 

-.0023 

.0042 

.0062 

.0090 

10. 

0. 

-1C. 

-15. 

-20. 

.05 

-15. 

7100. 

24. 

0. 

.645 

.0651 

.0C222 

-.5 

6. 

.8 

600. 

27. 

0. 

.1.5 

4.3 

-.0525 

1. 

10. 

l. 

0. 
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T'*LE  A-2 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


STRIP  ROTOR  THE  L'v  STABILITY  DERIVATIVES 
LIST  OF  INPUT  CARDS  {AFTER  130  AIRFOIL  SECTION  DATA  CARDSI 


CD 

COL 

SYMBOL 

1 

0-8 

TC 

COLLECTIVE  PITCH 

DEG  . 

1 

9-16 

OKU 

TIP  SPEED  RATIO 

1 

17-26 

VT 

TIP  SPEED 

F .P . S . 

1 

25-32 

S1G 

ROTDR  SOLIDITY 

1 

33-60 

AK 

ROTOR  ANGLE  OF  ATTACK 

DEG. 

1 

61-68 

03 

DELTA  3 

DEG. 

1 

69-56 

B 

NUMBER  OF  BLADES 

1 

57-66 

R 

ROTOR  RADIUS 

FT. 

1 

65-72 

RHO 

DENSITY 

2 

0-8 

OS 

SHAFT  ANGLE 

DEG. 

2 

9-16 

XR 

HORIZONTAL  ARM  C.G.  TO 

SHAFT 

FT  . ♦ 

2 

17-26 

HR 

VERTICAL  ARM  C.G.  TO  SHAFT 

FT  . - 

2 

25-32 

ET 

TAIL  LENGTH 

FT  . 

2 

33-60 

VTT 

TAIL  ROTOR  TIP  SPFED 

F.P.S 

2 

61-68 

S1GT 

TAIL  ROTOR  SOLIDITY 

2 

69-56 

RT 

TAIL  ROTOR  RADIUS 

FT. 

2 

57-66 

SVI 

AREA  VERTICAL  TAIL 

SQ  .FT. 

2 

65-72 

SHT 

AREA  HORIZONTAL  TAIL 

SO  .FT. 

3 

0-8 

VOLF 

FUSELAGE  VOLUME 

CU.FT. 

3 

9-16 

CDA 

FUSELAGE  DRAG  AREA 

SQ  .FT. 

3 

17-26 

EC 

ROTOR  BLADE  HINGE  ECCENTRICITY 

EC/R 

TYPE  OF 

ROTOR 

EC*0.  for  teetering 

UP 


EC*GT  . C.LT  1 . FOR  ARTICULATED 
EC«2.  FDR  RIGID 

EC-3.0  FOR  RESTRICTED  TEE  TER  1 NG 
6 25-32  TORSI  ON  SPRING  RATE  FOR  EC*3.0 

D3  IS  DELTA  3 FUTURE  USE  NON  0 
D1MENSI0NC0R(  18,5I,CDTU8,5) 

D I RE  NS!  UNCI.  1 (91 , 5 >,CD  1 19  1,  5 ) 

DIMENSION  Cl  (2  I. C2  (2  l,E  3 (2 > 

DIMENSION  E8I10I  ,29(101,(20(3) 

DIMENSION  CCU2,7),UX(12.7>  ,CBU2,7>  ,C  VT  (12  .7 1 ,CHT  (1  2,  71  ,CTR  ( 12 .7  I 
DIMENSION  CIO (2)  ,Cll(2»  , 02  (2  ) ,C  13  (2  >,CU(2  l.U  5 (2  ), Cl  6(2  1,(1 7(2  ) 
D013J«1.5 

READ  (5*1  ) (CL  1 (I  ,J  1.1*1,9  1) 

REA0(5«l J (CD  1 (1,J),I*1,9  1) 

READ(5.1)iCDT(l,J),I*l,lB> 

READ( 5 .1 1 (CDR ( I, J), 1 = 1,16) 

WR1TE(6,11(CLKI,J),1*1,91  ) 

WRITE (6*1  ) (CD1  <l,J),I*l,91) 

WRITE  (8,1  ) (COT  (1,J),I*1,  18) 

WR!TE(6,1)(CDR(I,J),!*1,18) 

13  CONTINUE 
WR I TE ( 6 ,6  > 

6 FORMAT (15Xt35HALL  COEFFICIENTS  BASED  ON  PlSC  AREA) 

WRITE (6, 5 ) 

20  RE  AD (5 ,2  ) TC ,OMU ,VT ,SI G , AK ,D 3 ,B , R, RHO 
READ(5»6)DS,XR,HR»ET,VTT  »SIGT  »RT ,SVT  ,SHT 
R E AD( 5 ,7 ) VOLF  .CDA.EC.RK 
9 FORMAT<6F8.2,F8.Q*F3.3,3F8.1) 

7 F0RWAT(3F8.2,E15.8) 
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WRITF  (6*5  ) 

WRITE  (6,2 )TC.0MU,VT,SIG,AK,D3,B,R,RH0 
WRITE (6,4)DS.XR,HA,ET,VT1,SIGT.RT,SVT,SHT 
WR I TE (6 ,7  JV3LF .CDA.EC ,RK 
WR  I! E (6  ,1 52) 

DS*DS/57  *3 
V*OHU*VT 

AS=-12.  $DAS*=4.  ISK*SK/57 .3 

S V * — l 5 . ADSY*15  . ATB*TC/57. 3 
J3*  7 

1 FORMAT (9F8»5 I 

2 FORHAT<2F8.5,F8.0,2F8.5,2F8.1 ,F8  .2.F8.5) 

3 FORMAT (3F8.3.5F10. 5) 

5 FORMAT!/) 

RC*SIG*3.141*R*R/B 
B1N*0  . 

W8*2»28*RC**l,l64 

G*BJN/(RHO*3.141*2.7*¥T**2»R*«3) 

G1*-WB*B/(RHO*32.2*.50*3.  141*R*»3  ) 

CRX*RK/(RHQ*3.141*R*®3*VT*$2) 

G 2s  IRK 

IF( AKJ30.3l.30 

30  AS* AK-OAS  V J3*l 

31  CONTINUE 

C ALPHA  LOOP 

DOl 02  JN*1,J3 

AS* AS*DAS  AAA*AS/57.3 

CR*0.  ACP*3.  SCW*0.  ACU*0  ACV*0. 

Al *-l  « AB1*-1. 

SK  *3. 

DO  32  1*1,12 
DO  32  J — 1 . 7 
CTR  C I ,J)*0. 

C V T ( I ,J)*0. 

CHT (I , J)*0. 

CO ( 1 , J ) *0. 

32  CC(l.J)*0. 

DO  107  L 1 *1  ,2 
A1*A1M  . 

SM*Al/57.3 

CALL  NO(TB,OMU.VTfSIG,SK*SM,AA.CLI  .COR.CMYSS.CHXSS.CYSS.COSS.CTSS, 
ICHSS.COSS.CLSS.CR.CP.CW.CU.CV  .D3.6.CD1.CDT  , SY  , E C ,G  1 ,R , G2  .HR  ) 

E9(Li )*SV 

C 2 0 1 L 1 > *CKYSS  SC11 (LI )*CMXSS  SC12( LI 1*C  YSS 
C 16(L 1 1 * C DSS 

C13(L1)*CTSS  ACl4(ll)*CHSS  SC  1 5(  L 1 )*CQ  SS  SC17(L1)*CLSS 
107  CONTINJt 

CMYDA  =C20  (2  J-C 20  (1 1 S C T DA  = C 1 3 ( 2 )-C  13(1) 

CXDA*-(C16(21-C16!l)J 

CMXDA*C11(2)-C11 (1)  SCHDA=C14(2  )-C14(l I 
CYDA  = U2(2)~C12(  1 1 SCGDA=C15<2  J-C15(l  ) 

SM  = 3. 

1 F ( CM  XD  A . NE  . 0 . I SM*-C  11  ( 1)  / ( 5 7,  3*CMXDA  I 
DO  108  LI *1,2 
n i -e  i ♦ i . 

SK*fcl/57.2 

CALL  NO  < T 8 ,OMU  ,V  T ,S  iG,  5 K ,5M  , A A,  C 11  ,C  DR,  C M YS  S,  C M X SS  ,C  VS  5,  C DSS  , C ! S S , 
1CHSS,C0SS,CLSS,CR,CP,CW,CU.CV,D3 .G.CDi.CDT.S/.EC ,G1,R, G2 .HR ) 

C20(L1  I=CMYSS  AC  11  (LI  )*=CMXSS  AC12{IU  = CYSS  JClfc(L!)=CDSS 
C131L1  ) *C  TSS  AC  14  (LI  )*CHS5  SC  1 5(  L 1 )*CQ  SS  SC17(L1>*CISS 
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108  CONTINUE 

CMV  DB*C20 (?) -C20( 1 1 ICTD6*C13I2 )-C13tl) 

CXDB*-(C16(2)-C16(1  )) 

CMXDB^CU  821  — C 11(1)  ICHDB»CI4<21-C19(1) 

C YDB*C12 (2  J-C12  (1  ) SCQDB-C  15(2)  -C  15(1) 

SK*0  • 

IF ( CM YOB  »NE . 0 • ) $K*-C20(  1)/(57.3*CMYDB) 

CALL  NO (T  B.OMU  *VT  »S  IG.SK .SM.AA.C 11 ,CDR,C MYSS.CMXSS ,C VS S, CDSS.CT SS . 
1CHSS.C0SS  ,CLSS.CR,CP.CW,CU,CV,D3,G,C01,CDT,SV,EC ,G1,R,G2,HR ) 
C10(1)«CMYSS  SC11  (1 )*CMXSS  SC12 ( 1 )*C  VSS  SC16(1)*CD5S 

C 1 3 1 1 ) «CT55  $C19(1)*CHSS  $C15U)*CC$S  U17(11»CLS$ 

81 *SK*57.3  AA1«SN*57.3 

WRITE (6,166) 

166  F0RMATI5X  .29HFINAL  VALUES  IN  PREVIOUS  LINE) 

WR I T E( 6 *1 52) 

WRITE (6,150) 

150  FORMAT  (lX,5HS.(AFT,4X,4HCnLL,6X.3HL AT  .5X.5HP1TCH  ,4X  ,3H|  IP  ,9X,5FSPEE 
ID) 

WR 1 TE (6  «1 51 ) 

151  FORMAT  (IX  .5HANGLE ,4X, 5HP  ITC H,  5X,  AH TR  1M,9  X,9HTR l M ,9X. 5H SP FED ,4X,5HR 
1 AT  I □ I 

WRITE (6,10 )AS,TC,A1  ,Bl  ,VT,OHU 

10  FORMAT  (9F8.2,F8.0»F8.3) 

11  FORMAT  (15X»7Ej.5.6) 

WR I TE (6  ,1 52) 

152  FORMAT!/) 

15  FORMAT! 10X.13HIN  ROTOR  AXIS) 

WRI TE (6,1 5 ) SWR1TE  (6,152) 

WRITE (6,153) 

153  FORMAT (5X,7HVAW  MOM ,6X ,9HP I TC H MOM ,6X,8HR0LL  M0H,6X,10H  Y , 
18X.6H  l ,8X  ,5  H X , 1 OX  ,6H  TOR  OUE  I 

WR  I TE  ( 6 ,1 59) 

159  FORMAT ( 5 X ,29HLAT  CYCLIC  EFFECTIVENESS) 

WRITE (6,11  )CMYOA,CMXOA,CYOA,CTDA,CHDA,CQDA 
WRITE (6,155) 

155  FORMAT (5X,26HPITCH  CYCLIC  EFFECTIVENESS) 
WR1TE(6{11)CNYDB»CMXDB«CYDB»CTDB .CHDB.CODB 
CC ( 1 ,1 ) “CMYOA  S CC ( 1 ,2  ) =CMXOA  SCC ( I ,3 i = C YD A 
CC ( 1 ,9 ) «C  TOA  SCC ( 1 ,5 ) *CXDA  SCC (1  ,fc)*CODA 

CC (2,1 ) -CM YOB  SCC  (2,2  1 *CMXD6  SCC ( 2 ,3) r CYDB 
CC ( 2 ,9 ) *CTDB  SCC  (2.5  )«CXDB  S CC ( 2 ,6) *CQD B 
J1  =2 . 

TB1*TB*1  ./57.3 
D099J J*1 ,6 
Ji^JlU 

CR=0.  SCP*0.  SCW=0.  SCU=0.  $CV*0. 

IF(JJ.EO.l)  CR=l./57.3 
IF(JJ.E3.2)  CP=-l./57.3 
IF(JJ.EQ.3)CU*-1  . 

1 F ( JJ »E0 ,9 )CW*1 . 

IF ( JJ»E0.5)CV»*l . 

IF  (JJ.EL.6)TR*TB1 

CALL  N3 ( TB. 3MU, VT, S IG, 5K. SM, A A, C L1,C DR, C MYSS.CMXSS  ,CYS S, CD55 ,CT  SS , 
1CHSS«CQSS,CLSS,CR,CP,CW,CU,CV,D3,G,CD1,CDT,SV,EC  ,Gl,a,G2,HR) 

F B( JJ ) -5V 

CIO  C 2 ) = CMYSS  SC11(2)=CMXSS  SC12(2)=CYSS 
C16(2)=CDSS 

C13(2)=CTSS  SC19  12  1 *=  C H SS  SC15(2)=CQSS 

CMYD*C10  (2  >— C 1 0 ( 1 ) SCT0-C13C  2)-C13<  1) 

CXD=-!C16(2)-C16(1) ) A~ 


cmx o* ci i(2}-cii(i)  $cHD*tum-mm 
CYD«C12(2)-C12(1  ) SCQD  *C  1 5 t 2 )-Cl  5 ( 1 ) 

WRITE  (6,152) 

WRITE (b* 1 S3  ) 

IF(JJ.EQ.l)  WR IT  E (6  ,157 ) 

1FOJ.EQ.2)  WR  IT  Etb  ,158 ) 

1FIJJ.E0.3)  WR  IT  E(6  *160  ) 

1F(JJ.E0.4)  WRITE (6*159) 

IFUJ.EQ.5)  WRPE(6.161) 

1 F ( JJ  ,E C .6  )WR  1 TE  (6,163) 

JF(JJ.EQ,6)TB*Tbl-l  ,/57.3 

WRITE  (6,11  KMYP.CHXD.C  YD  .CXD,  COD 

CC  f Jl  ,1)*CMY0  SCC(  JI  ,2  ) *CMXD  SCC ( Jl ,3 ) *C YD 
CC(  J1  ,4)  *CTD  SCC(J1,5)<HD  ICC  ( Jl»6 J -CQD 
IF (JJ.E0.6J60T099 
!F(JJ -4)99, 162, 162 

162  CMY0*CMYD«V  SCMXD*CMXO*V  SCYD*CYD*V  SCTD*CTD«V 
CHD*CHD«V  ACQD=CQD*V 

CXD^CXD^V 

lF(JJ«Efa,4) WRITE (6,164  ) 

1 F ( J J . EQ . S )WR  I TE  (6.165) 

WRITE  (6,11  )CMYO,CNXD,CYD,CTD,CXD,CQD 
J1«J1*1 

CC  ( J1  ,l)*CHfO  SCC(  Jl,2)  *CMXD  S CC  ( J 1 ,3  ) =C  YD 
CC  ( J1  ,4 ) *CTD  SCC(  J1,S)«CXD  SCC  t Jl,6)  -CQD 
99  CONTINUE 

163  FDRHAT(5X,10HC0LLECTIVE  ) 

164  FORMAT  < 5 X .1SH4NGLE  OF  ATTACK) 

165  FORMAT (5X  ,8HS  IDESLIP) 

157  Fl)RMAT(5X  ,12  HR  OIL  DAMPING) 

158  F0RNAT(5X,13HP1TCH  DAMPING) 

159  F0RMAT(5X,17HVERT1CAL  VELOCITY) 

160  FORMAT (5 X , 1 6HF OR  WARD  VELOCITY) 

161  FORMAT (5X.13H5IDE  VELOCITY) 

WR 1 T E ( 6 ,1 52)  SWR  IT  E(6 ,40)  S WR 1 TE (6  ,152  I 
WRITE (6,150) 

WR 1 TE (6,1  SI  > 

WRITE (6,10  )AS,TC, A 1,81, VT.OMU 
WRITE  (6,152) 

WR1TE(6,153) 

0041  1*1,4 
D041  J*l  ,7 

41  CCt 1,J)*CC(1»J)»57.3 
D042  J *1 ,7 

42  CC(10,J)*CC(10,J  )*57.3 
CALL  PRNT  (CC) 

40  FORMAT (10X .15HSUMMARY  UF  DATA) 

A54*CC(  5 ,4)  A A64*C  C( b ,4  ) SA55*CC(5,5)  SA65=CC(6 

CHANGES  ROTOR  DERIVATIVES  TO  STABILITY  AXIS 
CC(5»4) =-(A54-A64'DS)-(A55-A65»DS)*DS 
CC(5,5)=-(A54-A64*DS)«DS'*(A55-A6  5«DS  ) 
CC(6,4)=-(A64*A54"DS)“(A65*A55«DS)$DS 
C C ( 6 ,5  ) * - ( A64*  A54*DS)*DS-*(  A65*A5  5*DS  ) 

CC(5,1)  = -HR*CC(5,5)/R-XR<*CC(5,4)/R 
CC(6,1)=-HR*CC (6,5)/R-XR»CC (6  ,4 > /R 
CC (1,4  ) =-CC (1,4) 

CC(2,4)=-CC (2,4) 

CC ( 3,4) =-CC (3*4) 

CC(4  ,4) = — CCt  4 ,4) 

CCt  7 ,4)*— CC(  7,4) 
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CCI8,4)«-CC|8,4) 

CC  (9,4 ) *-*CC  (9,4) 

CCU0,4)*-CC(10,4) 

CCI  1*21* CC(  1 ,2 )*CC(  1 * 31  *HR/R 
CC  C 2 *2  1 *CC(  2 *2  ) ♦ C C C 2,3)*HR/R 
CCI  3*21 *C C C 3.2MCCI  3,3)*HR/R 
CCI  9 *2  ) *CC  ( 4.2KCC!  4,3)*HR/R 
CC ( 5 *2  ) *CC { 5*2  1 ♦CC ( 5,3)*HR/R 
CCI  6*21  *CC ( 6.2MCCI  6,3»*HR/R 
CCI  8 ,2 ) * CCI  8 *2  5 ♦ C C I 8,3)*HR/R 
CCI  10  *2  ) *CCt  10  ,2)*CC(10,3)»HR/R 
CC(11«2)*CC(11«2HCCI11»3)  *HR/R 
CCI  1,1)*CC<  1*1)«CCI  1 »4)*XR/R-CC  I 1*51  *HR/R 
CCI  2*1)  *CC  I 2*1  )*CC  I 2»4)*XR  /R-CC  I 2.5MHR/R 
CCI  3.1 1 *CC  I 3,1)  ♦C  Cl  3,4)»XR/R-CCI  3.5MHR/R 
CCI  4 *1 ) *CCI  4 *1 ) ♦ C C I 4,41  *XR/R-CC  I 4,5)*HR/R 
CCI  5 *1  )*CCI  5,11-KCI  5,4>«XR/R-CC I 5,5)*HR/R 
CCI  b « 1 1 *CC  I b.lMCCI  6 ,4 ) *XR  /R-CC  ( 6,5)«HR/R 
CCI  8,1  1 *CC  I 8*1  MCC  I 8,4)  *XR  /R-CC  I 8*51  *MR/R 
CCI  10.1  )*CC  110,1  )*CC(10.4)*XA/R-CC(10,5)  *HR/R 
CCI  11  ,11*CCI11  ,1 1 •* C C 1 1 1 ,4) *XR /R-CC  I 1 1*5 ) *HR/R 
DO  14  1*1  ,10 
14  CCI1, 71*0. 

CCI10, 7) *-CC <10,3  )*XR/R 
CCI1,7)«-CCU,3)*XR/R 
CC  (2»7)*-CCI2,3)*XR/R 
CHANGES  10  COHPIETE  NOOEl 
TAIL  ROTOR 
OMUT*ONU*VT/VTT 

16*0.  $SK  *0 . 1 SM*0.  I AA*0.  $CR*0.  $CP«0.  SCW»0. 

(U*0.  SCV«0.  $03*0.  SEC *2. 

CALL  NOIlB.OHUl.Vll ,SIG1 ,SK ,SM,A A, CL  1 ,C DR.CHY SS ,C*XSS, CV SS.COSS »CT 
1 SS, CHSS, CQSS, C LSS, CR, C P, CW, CU. C V. 0 3, G, CD  1, CDT.SV.eC, Cl  ,R  T ,02 ,HR  1 
CT6*CTSS 
TB*l./57.3 

CALL  NO (TB.QMUT ,VTT,S16T ,SK *SK,A A,(L l.CDR.CMVSS ,CMXS S, CV SS ,CDSS , Cl 
1 SS,CHSS,COSS,CLSS,CR,CP,Cil,CU,CV  .0  3,0, CD  I, CD!  ,S  V ,EC,  01  ,RT  ,02  , FR  1 
C1IB*(CTSS—CTB 1*57.3 
CW*1  . 

TB*0. 

CALL  N0ITB,0HUT,VT1,5ICT,SK,SM.AA,CL1,CDR,CMVSS,CMXSS, CVSS.CDSS  ,CT 
1SS«CHSS,CQSS,CLSS«CR,CP,CV,CU,CV  ,D  3,O,C0  l.CDT  ,SV  ,EC,01  *RT,G2,HI) 
CTW*C1SS-CTB 
TRR*(RT*VTT/IR  *VT 1 1 **2 . 

YVI  *-CI  W *TRR  $ZVT  *-Y  VT*  El/  R SYBTT*  V*VVT  SZBT1*ZVT*V 

YTCT*CTTb*TRR  SZ1 C 1 *-Y TC T*E 1/R 

YSDT*-ET*YVT  SZ  SDT  *-Y  SD T*E  T/R 

HT*5. 

YTD! *YVT *H1 
ZTDI*-Y1  DT  *ET/R 
C VERUCAL  TAIL 

CLAT*2.0  $8C=VT»VI«3.141*R*R 

E VT  *CLAT* V*V*S VT /2 . SYBTVT *-EVT /QC 
ZBT VT  *0.  $Y VV1  *0.  $ZWT*0. 

IF  IV.LT.l . 1 G0T08 

ZETVT*£VT *EI/IOC*R>  $ YV VT =YB 7VT /V  SZVVT'ZBTV ' /V 
8 CONTINJE 

YSCVT*-YVVT*ET  IZSVDT =-ZVVT*ET 
C HORIZONTAL  TAR 

CLAT=2.t  a-23 
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EHT »C  LAT  • V*V*$HT/2. 

ZAHT«-EHT*(1  .-(E8(9)-E9(im/QC 
2WHT*0. 

2TH0T  *0 • 

If  (V.GT.l  . I ZWHT *ZAHT /V 
IF  < V.GT . I . I ZTHDT *-EHT  *ET/ (QC  *V 1 
CMAHT«ZAHT»ET/R 
CMWHT  *ZWHT  *ET/R 
CMTDHT=ZTHDT*ET/R 
C BUOY 

CXU*-CDA*V/(  R >»VI <*VT  ) 

HE  *.8 

CHAF«2.*HF*V0LF*(1  .-(E8(9)-E9I1)  1 » /(3.19  1*R«*3  . I 
CNBF*2 .*HF*VOl F/f 3 .191 «R*« 3 .1 
CHAF«CMAF*V»V/(2.*VT*VT) 

CNBF«CNBF*V*V/(2.*VT*VT  ) 

CMWF  *0. 

CNVF  *0. 

IFU.GT.O.)  CHW  F*CH  AF/V 
I F < V • OT .0*1  CNVF»-CNBF/V 
CTR ( 8 ,3  1 *YVT i CTR ( 8 *7 1 *Z VT 
CTRI9.3  l*Y8f T*  CTR 19,7  >«2BTT 
CTR (12*3  )*YTCT$  CTR  <12,7  1«2 TC T 
CTR (l 1 »3 ) «Y$OT I CTR  «11 . 71-ZSDT 
CTR(3»71  »YVT  *HT 
CTR  C3  * 7) *ZTDT 

CTR « r,,2J«-'TRf  3,7)HT/n 
Cist  *;l  j«-CTR  ( 8,7  )<MT/FT 
CTR I 9,2  »«-CTR  I 9,7}*HT/ET 
CTR ( 1 1 ,2 ) «-C  TR (ll,7)*HT/ET 
CTR(12,2)*-CTR(12.7)*HT/ET 
CVT(8,3)«YVVTS  CVT (8,71*  ZVV1 
CVT (9,3) *YBT  VT S C VT  (9. 7 1 *ZB7VT 
CVT (11,3 )«YSDVTS  CVT (11, 71*ZSVDT 
CHT (7,9  )«ZAMTS  CHT (6,9 1 = ZWHT 
CHI (7,1) *CHAMT  $ CHT (6 , 1 1 *CMWHT 
CHT  (9,11  *CFTOHT  I CHT(9  ,9  l*ZTHDT 
C B( 5 ,5 ) *C  XU 

C VT(  8,21 =CVT ( 8 ,7 1 *HI / ( 2. *ET  1 
C VT  ( 9,2  1 *~C VT  ( 9,7  )*HT / (2,*E  T 1 
CVT  (11,2  I— CVT  (11,7  l*r|T  / (>  *#E  T 1 
CVT  ( 3,31  *CVT  ( 8,3)*HT/2. 

CVT ( 3 ,2 1 *CV T ( 3,3 >*HT/(2«»R) 

CVT  ( 3,7)  *CVT  ( 3,3l*U/(-M 
CB  ( 7 , 1 1 sCMAFS  CB  (6.11-CMW 
CB  (9,7  l^CNBFS  CB  (8 , 7KNV‘ 

WRITE  (fc, 152  IS  WRITE  (6,  IMS  WRITE  <6,  152) 

WRITE  (6.1531S  WR  1TE(6,152) 

33  FURMAT(5X .18HMA1N  ROTOR  #>E R DEG) 

39  FORMA? (5X ,9H800Y 1 

35  FORMAT  (5X.13HVERT  1 C A>  TAIL  1 

36  FORMAT (5 X , 15HH0R l ZONf AL  TAIL) 

37  FORMAT (5X .10HTAIL  ROT 3R 1 

38  FORMAT (5X ,16HTCTAL  PER  P AD  I AN  1 
WR  I T E ( 6 ,3  3 1 
CALL  PRNT(CC) 

WR 1 TE  (b  ,39  ) 

CALL  PR \ T (C  3 1 
WRITE (6,35  ) 

CALL  PR  M (CVT) 
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UR  1 IE (6*36) 

CALL  PRNT(CNT) 

WRITE  (6,3?) 

CALL  PRNT(CTR) 

00  39  1*1,12 

00  39  J*1  ,7 

39  CC<  1,J)*CC(  I , J)  ♦ce  I 1,J)*CVT(I ,J)  ♦CHTII.J  )+ClR  (J  , J> 

W ft  1 TE  (6,381 
CALL  PRNT  (CC  ) 

102  CONTINUE 

16  FORMAT ( 10X  ,1 7H IN  STAB1L  1 *Y  AXIS) 

C0T020  SEND 

SUBROUTINE  N0(  T8,0MU,  VT , SI  G , SK,  SH,  AA  ,CL  1 ,CDR,CM V SS  ,CHX  SS  ,CYSS  »CDSS 
1 ,CTSS,CHSS,CQSS,CLSS,CR ,CP ,CH  *CU ,C V, 03,0 ,C0 1,C0T  ,SV , EC ,G 1 ,R ,02 ,HR ) 

0 1 MEN  SI  ON  COR (18, S), COT (18, S) 

DIMENSION  CL  1(91,5),  CO  1(91,  5) 

DIMENSION  C1(2),C2(2),E3(2) 

DIMENSION  C10(2),C11(2),C12(?),C13(2),C1<>(2),C15(2),(16(?),C17(2) 
DIMENSION  86(24), B7(24), 810(24) 

DO  8 J7 *1  ,24 
96 ( J 7 )*0. 

B10(J7)*0. 

8 B7(J7)«0. 

CRK*G2 

J*1 

DT*VT*57.3/(R*15.  1 
0T*1./0T 
152  FORMAT!/) 

DSY*1 5 • $SY«-15. 

T C*TB*57 .3 
AS>AA*S7.3 
C *1 120. 

SMF-.014 
S V =5 . 

8LAM-0. 

S F*OMU*S INI AA) 

S2*  SF  *SF 

02  *OMU*QMU 

INDUCED  VELOCITY  LOOP 
4 FORHAT(3F10.3,4F10.5,F10.3) 

DO  101 JJ*1 ,5 

1 F( SV .GT • 1 13 . ) 5V=100.«.1»SV 
D099I J*1 ,2 
E 1 =SV/VT 
E2  *E1 *E 1 
E3(1J)*E1 

C2(U)  = SQRT(4.*E2»(02*S2-2.*SF«E1*E2  )) 

IFfSV.LT.O.)  C2( 1 J)*-C2 ( 1J ) 

14  F0RMAT(4X,2HAS,7X,2HSY,8X,3HAl2teX,3HCLU  .7X.3HCDU.6X  *4HSPAN  ,7X  ,3HS 
1TM  ) 

IFIJJ-5 >22,21, 22 

21  WRITE (6,14) 

22  CONTINUE 
IFU.EO.DGU  TO  it 
DO  10  J7-1  ,24 

N1  =J7-1 

IF(Nl.EU.O)  V 1 =24 
N2= J7+1 

!F(N2.EL.25)N2=1 

B8  = (B7(  J7 ) — 3 7 ( N 1 ) )/UT  A'25 


B9«(B7<N2  )-B7(J7  ))/0T 

10  Bfc  ( J7  )* (B8*B9  ) / 2 • 

00  11  J?*l,24 

N 1* J7-1 

IF(Hl.EQ.O)  HI  *24 
N2*J7M 

IF (K2.EC.25)N2*1 
B 8 * (Bb  ( J?  )-B8  (N 1 1 )/DT 
B9*(8b(N2)-,B6CJ7) I / D1 

11  B10I J7)«<88*89)/8. 

12  CONT’fJ'JE 

0 01071**1 ,2 
C1SS*0.  ACHSS*0. 

SY—18. 

C AZIMUTH  LOO? 

CKXSS  *0  • ACHY  5S*0  , ACYSS*0, 

CTSS^O.  ACDSS*0.  acoss*o. 

D01  00KJ*1  ,2!i 
S Y*SY  ♦DSY  $SA*SY/S7.3 
DTM*0.  1 B 1*0. 

BQ*0  • 

AR*  EC 

1 F(  EC. GT . 1 . 1 A R*0 . 

CTS»0.  ACPS*0.  »CHS*0. 

CHX<*0  • ACMY*0.  ACYS*0. 

CRS  «0 . 

C SPAN  LOOP 

DTH*B7<  JK)*D3/(87.3«8.I 
BP*57.3*87(KJ) 

IHU03K<  *1 .8 
DAR*U.-EC  1/8. 

IF  CEC.5T.1.)  0 Ak * . 1 28 
lF(RK.FU.l)  DAR=DAk/2. 

AR  *AR  + DAR 

I F( AR . F 0. 0. I A R* ,0625 
C FOR  THO/OHEGA  RA1ES 

VY*0HU*V7»SIN(AA  )-SVm«R  • (C  R*S  IN  C S A > ♦ C P *CDS  i S A ) ) HW 
VY*VY*AR*R«B6t JK  ) 

VX*VT*(  AR*0MlJ*SIN(  5AMC0SI  AAJ  }«CU«SIM5A  >*C  V»C 3 8 ISA  ) 

STM* ABS  C V X > / C AIM*  STM/ , 2 
1 F ( 1M)2B.27,28 

27  1 H *■!  I PM  *0.  A 1 MM  = 2 

G0T029 

28  IMH*IM*1  AAUM*1N  $AJM=Al»M«.2  APM* (STM -AUK) / .2 

29  CONT  IKUE 

G A* ABS ( V Y /VX ) 

C AH  LOCAL  ANGLF  OF  ATTACK  RELATIVE  WIND  TO  PLANE  PE  RP IN  D1 CUL AR  TO  SHAFT 
AH* AT  AN (G A ) 

IF(VY)41,40,40 

90  I F(  VX  ) AN  ,93,93 

44  AH=3.  141-AH  AG0T043 

91  IF  (VX  196,99,48 
48  AH=-AH  t GOT04  3 

46  AH* AH-3 .141 
43  AH*  AH 

T T = -8  . ~ AR  / 67 . 3 $FN=1. 

SB  =SA-D  3/87. 3 

All=TE*6./6  7.3*TT-SM*S!N  (SB  )•»  SK  * CO  S ( SB  I - CTH-*AH 
IF  UL1  )4?  ,47,49 

47  A L 1 *- AL 1 iFN* -1  . 
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49  AL2*AL1»57.3 
ID*AL2/2.*1. 

AUS*I  0*2 

PIW.-IAU5-AL2  )/2. 

C CLU,  C0J  LOCAL  LIFT  AH  DRAG  COEFFICIENTS 
1 D0*ID*1 

CL3«{  CL1I  10,IM)«Hi«CCLll  iDD , I K l-CL 1( ID, I H ) ) >*F N 
CD3*l CDI I JO,  IMMPU*(CD1(10D»U  l-COl  1 10,1*1)1 
CL4«FN*  (CL  HID,  I MM  )*PU*<CL  1 (I  UD,  1HM  1-CLU  ID.IHH  II I 
CD4-CD1  UO.IHH  l+PU*  (C3  1 1 IDD  ,1  MK  ) -CD1  UD.IHH)) 

IF  C 10-171 33,33 ,37 

33  IFIKK -6)34,34, 35 

34  PTC-l-KK/6 

CD5*CDR(I0»IM)*PU*CCDRIIDD»IH  )~C0R(10,»H  )) 

CD6<DR  (1  D.IHH  )*PU*  (COR  1 100 ,1 NM  ) -CDR  (IO,!MH)) 

GO? 036 

35  P7C*(K5-6)/2 

C05>CDT( IO»IH)«PU*((DTI I DD , 1H )-C 01 ( 10, IH  )) 

CD6»CDT  ( I 0»IMH )•♦  PU*  (CD?(  ID D * 1 HM ) — C D T (ID,  IMM  ) ) 

36  CD3*CD3*PTC»(Cl)5-CD3) 

CD4*C04*FTC*  IC06-C03) 

37  CONTINUE 

CLU*CL3*PH«ICL4-CL3  ) $CDU»CD3*PM*  IC04-CD3) 

C CX,  CZ  COEFF  PARRALEL  AND  PERPINDIC  UlAR  TO  BLADE  TWIST  AXIS  IN  PLANE 
C PERPINOICULAR  TO  SHAFT  AXIS 

CX*(-CLU*S1N  IAH  )4C0U*C0S  IAH  )1 » I VX» VX  *VY*VY ) *5  I C / VT 
C2*(CLU»C0SIAH )4C0y*SlN(AH ) )• (V < *VX* VV «V Y) *SI &/ V T**2 
EA*EC 

JFIEC.&T  .1.1  EA-AR 
1FIEC.E0.0.)  EA*HR/R 
DC  Z *0 . 

DCX*0. 

S9*SB 

CDSY*-CX*C0SI5BI-tZ*S!N(B7{KJ,-)‘  SI  N ( S9) » OCX «C3S  ( SB  > 
C0SC‘CX*SINISb)-CZ*SIN(87(KJ)  )*l  OS  C S9)-D  CX*  SI  N ( SB) 

IFIEC.EO.O.)  GOT 023 
G0TD24 

23  EA=EA*S!N  (B7IKJ)  ) 

IF  (KJ.GT.12.  ) EA*-EA 

24  CONTINUE 

CHXC«-CZ*C0:(8  7IKJH«EA*SIN(SEI 
IFIFC.EQ.O,)  GOT 025 
GOT  026 

25  EA«EA»SIN  (B7  (KJ)  } 

IF IKJ.GT.7.0R.KJ.L!  .19)  EA =-E A 

26  CONTINUE 

CMY  C*-C2  *COS ( B7I KJ) )*E  A® CCS! SB) 

CYS»CYS*CDSY 

CMX*CHXHHXC  KHY*CHY«CMfC 
C Z =C  Z-DCZ 

CTS<TS*CZ  KH5*CHS*C0SC  $C  PS=CPS*CX  *AR 
CRH*-CZ*C0S(B7CK J) I *E  A 
CRSsCRS*CRM 
IF(KK -9)106*104, 106 
104  CTS=CTS-CZ/2. 

C Y5=C YS-C  OSY/2. 

CHS=CHS-C0SC/2.  $C»-S-CP:-CX/2  . 

CMX=CHX-CMXC/2 . *C HY  = C HY-CHY C/2 . 

CRS-CRS-CRM/2 . 

106  CONTINUE  A'i>7 


A S =AA*57 . 3 
AL 2 ==FN*AL  2 

Bl=Bl4C7/(Gl»ARl-32  .?»R/ (AR*VT*VT1 
lF{£C.e«.2ai  61=0. 

I Ft JJ-5 1 1 9,5 ,1 9 
b CONTINUE 

IF (LK-2  119,20, 19 
2C  CONTINUE 

lFUJ.EO.l)  &□  TO  1 7 
lFiKJ-7118,17,18 

18  IFCKJ-191 19,17,19 

17  WR1 TE (5,4 >AS,$Y,At2,ClU.C0U,AR  , SIN, BP 

19  CONTINUE 
103  CONTINUE 

CT *CT S/8  . $CH«CHS/8. 

C7*(HX/8 . $C8=C*Y/8. 

CHX0S«C7/2.  »<  HYOS«C8/2. 

CTOS-CT/2. 

CDDS-CH/2. 

CY0S*CYS/16. 

CRY  *CRS/ 1 6 . 

IFIEC.EQ.3.1  bl=CRY/CRK 
CQOSKPS/16. 

B P*81  *57. 3/8  . 

BU*bl  /8 . 

B7(KJ)»BC 
B 1 *0 . 

SJ«SK*57.3 

15  FQRNAT(4X,2HAS,8X»2HTC»8X,2HTIl»7X»4rtCT0S  ,6X  ,4HC  DOS  ,6X,  AH  COOS  • 7X  , 2H 
1SV) 

105  C T SS*C  T SS  «C 1 3S  *C0 SS =C 0 SS *COOS  SCQSS»CQSS*C03S 
CMXSS«CMXSS*CMXG5 
CHVSS«CMVSS*CMYUS 
CYSS*CYSS*CV  OS 
100  CONTINUE 

IF  IEC.EQ.3.)  GOT 07 
IF  (EC  .GT.0. 1 GOT 05 
7 CONTINUE 
DO  3 K J= 1 ,12 
KL*K  JM  2 
A8  * 1 . 

8 7 (KJ  1=A8  * (8  7 (K  J 1-8  7 (XL  1 1/2. 

3 B 7 (XL  1=-B  7 (X  J ! 

5 CONTINUE 
107  CONTINUE 

CTSS*CTSS/24.  $CDSS*:D  SS/24.  V CQSS  'COS  S/24. 

CMXSS=CMXSS/24.  $C  MYS  S *(H  Y SS / 25  • 8CYSS -CYSS/24  . 

Cl SS=CT55*C3S (AA 1-CDS5 *S1  N (A  A 1 
CHSS*CTSS«SIN(AA  >*C  DSS«-C  US  ( AA  1 
C TC  *C2 ( 1 J 1 
SS=SM*57.3 
FNC=G(4J*C1  SS/COSS 
WRITE  (6,152) 

WRITE  (6,16) 

WRITE (5,2 1 

WRITE  (6,9  >AS,TC,SJ.VT,UMU,CTSS,C  OSS,  CtSS.CrtSS.CQSS.SY,  CTC.ENC 
WR  1TE(6,1  1SS,CMXSS.CHYSS.CYSS,BQ,S8, AH 

16  FORMAT(4X,2H4S,5X»2HTC»7X*2HT8,bX»2HVT,4X»3H3l(J,5X  ,4HC  TS  S , 8X  ,4HC  0 S 
1 S,ex,4HU  SS.7X.4HCHSS,  7X  .4HCQS5,  7X  ,2HSV  ,7X  ,3HCT  C ,6X,3HENC1 

9 FORMAT  (3F6.2.F8. 0,F6.3,5F11  .6, FI  1 .3  , F 1 1 . 6,F  8 . 3 1 ‘V 


1 FORMAT  { F 6 . 3,  (»F  11 .6) 

2 FORMAT (2X,2HSS,5X,5HCMXSSt7X,5HCMYSS.7X»*HCYSS) 
Cl  4 1 J )*CTSS 

J9*JJ-5 

IF  (J9. EC. 0)6010102 
SMF  *SHF-.001 
99  5V*SV*$MF*VT 

SL-CCitl)  -C2  (1)-Cl(2)*t2l2))/IE3(1  )-E3(2  I) 

El«  ( Cl  ( 2 ) -C2  ( 2 I)  / SI 

101  SVME3(2)~El)*VT 

102  WRITE(6#152) 

RETURN 

END 

SUBROUTINE  PRNT(UX) 

DIMENSION  UX C 1 2 • 7 1 
W R I TE (6*21 ) 

1 F ORMAT (5Xs2AHl AT  CYCLIC  EFFECTIVENESS! 

2 FORMAT <5X,26HP1TCH  CYCLIC  EFFECTIVENESS) 

3 FORMAT (SX  »12HR0LL  DAMPING) 

4 F0RMATI5X  .13HPITCH  DAMPING) 

5 FORMAT ( 5X  1 16HF  QRWA RD  VELOCITY) 

6 F0RMATI5X.17HVERT1CAL  VELOCITY) 

7 FORMAT (5X , 15H ANGLE  OF  ATTACK) 

A FORMAT (SX.13HSIDE  VELOCITY) 

9  FORMAT (SX ,8MS1 OESL IP) 

10  FORMAT ( 5X ,16HR0T0R  COLLECTIVE) 

11  FORMATISX.llHYAW  DAMPING ) 

12  F ORMAT ( 5 X t21HT A 1 1 ROTOR  COLLECTIVE) 

17  FORMAT (BE  IS. 6 ) 

DO  20  1*1,12 
1FU.E0.1)  WRlTE(b.l) 

I F(  1 . EQ  .2  ) WRITE(6.21 
IF  I 1. EQ. 3 ) MR! TBlb.3) 

IF  Cl . £0.4  ) WRITE  (6. A) 

IF  (1.EQ.5  ) WRITE (6«  5) 

IFU.E0.6)  WR  I TE  (6*6) 

I F ( 1.EQ.7)  WRITEC6.7) 

I F ( I.EO.A)  WRITE(6.8) 

IF( I • EQ.9 ) WRITE (6,9) 

IF (I.EQ.IO)WRITE (6,10) 

IF  ( I . E 0. 1 1 )WR  1 TE  (6*  1 1 ) 

1F(I.EQ.121WRITE(6.12) 

20  WR1TE(6,17)(UXI1,7),IUX( I,J),J*1,6)1 

21  FORMAT!/) 

WRITE (6,21 ) 
return 

END 

* L I ST ( STOP  ) 
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APPENDIX  B 

PROGRAM  FOR  SOLUTION  OF 
LONGITUDINAL  EQUATIONS 


BASE  VALUES 


TABLE  B.l 

EXPANSION  OF  LONGITUDINAL  DETERMINANT 
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TABLE  B.l  (Continued) 
EXPANSION  OF  LONGITUDINAL  DETERMINANT 
PITCH  CONTROLLER  TERMS 


Note  triangular  relationship  between 


de 

d3l 

cT 


9 § 0 SAS  Contribution 


Ik 


B-3 


SAS  CONTRIBUTIONS 


TABLE  B.l 


EXPANSION  OF  LONGITUDINAL  DETERMINANT 
(Continued) 
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V 1 lT 


-C  Iy'C_  . C C c7  -C_  C.C 
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v ra0  ZB1  % Z§  mBl  mu  L ZB1  dBl 


iii-C  C 
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* Note  triangular  relationship  between 
a & SAS  Contribution 


^Quintic 


nnnon  oooooo  ooooooooooo 


i i t * ' * 


LIST (STA«T ) 


SOLUTION  Of  i ONfS Itij^ImAL  EOU A t I nMS 


THIS  PROGRAM  EVA.UATES  THE  LONGITUDINAL  STABILITY  OUARTJC 
COEFF I C I ENTs  A i B * C i D < p OBTAINS  ThE  RESULTING  ROOTS, PERIOD 

and  time  to  damp  to  half • amplitude,  to  obtain  satisfactory 

STABILITY  AUTOPILOT  GaInS  c50aBl/TDD,C5U8l/TD.  AND  CSOsBl/THPTA 
MAY  BE  US£D 

THE  REQUIRED  INPjT  CARDS  AnD  PROPER  COLUMNS  ALONG  WITH  THE  SympoLS 

_ - _ _ — . r-\  _ _ f~l  _ — _ . . _ _ r*  . . 


AND 

CARD 

1 


IDENTIFICATION  aRR  PRESENTED  BELOW. 


COLUMN 
1*15 
16-30 
31-45 
46-60 
61-75 
1-15 
16-30 
31-45 
46-60 
61-75 
1-15 
16-3Q 
31-45 
61-75 
1-15 
16-30 
31-45 
46-60 
61-75 
1-15 
1-15 
1 -30 
31-45 
46-60 
460-75 
1-10 


SYMBOL 

PM 

V 

CL1 

ALPH 

Y 1 
CXU 
CXA 
CZU 
CZA 
CZAD 
CZTD 
CMU 
CM  A 1 
CMTD1 
CZTC 
CMTC 
CZRl 
CMP1 
CM! 

CXI 

Fm 

GV 

C50 

C51 

C52 

MORE=l 


mass  m#v/c*ad 
velocity 

lift  COPFFICIPNT 

ANGuE  Op  ATTACK 

moment  nr  Inertia 
Stability  derivative 


FT/SbC 

W/Q*AD 


AD  DISC  AREA 


I yy/q#ad«Rotr 


DCX/DU 

dcx/dalpha 


tor  STABILITY  DERIVATIVES  see  HELIo 


COI LECTlv:  PITCH 

collective  pitch 
cyclic  PITCH 
CYCLIC  PITCH 

horizontal  tail 


autopilot  SaIn 
autdptlot  gain 

AuTOPfLOT  COEFFICIENT  S2 


control  effect 


gains,  -1  NEW  SFT, 


typical  InPuTS  aRf 
. 006? 


C, 000194  1.  .006?  o. 

C*  * 0000042  . 0000079  -’,  000009  -.000043 

CO,  ,000000078  .000000008  0. 

C--.041  .00000336  0.  -.0013 

C , 0043 

CO,  0,  0*  0.  0, 

Cl 

cl.  it  It  2.  4» 

c-l 

Cl, 26  100.  .405  0-  , 

C.  000486  . 05207  -.  000173  -.793 

C0 , .00035  -;o4?7  -.000 

C-3.14  .155  .793  --1°8 

C.353 

CO,  0,  0,  o’.  0. 

CO 

C PROS 

DIMENSICNxCOFdi)  ,COF(li ) , ROO^R  < io  > .ROOT  I < 10 ) 


. 00051 
0. 

.0000002? 

0. 


0 • 

- . 793 
- .00088 
-.108 


.033 

0. 

-.13 

-.0144 


COMMONXCOF.COFiM*ROnTR#ROOTI . TER 
M»4 

98  CONTINUb 

READ<5,6)PM,V.CL1.AlPh.Y! 

READ<5.6)CXJ,CXA,CZy,CZA»CZAD 
READ{5»6)CiTDl#CMu*CMAi.CMAn#rMTDl 
READ<5i6)CZTC,CMTC,CZBl.CMBl.PMl 
READ  ( 5 » 6 > CX  I 

1 FORMAT ( 7F10 . 6 ) 

WRJ1E(6»2) 

2 F0RMAT(t>X,48HAlRCRAFT  QA  INPUTS  IN  ORDPR  *Rb  PM  , V , CL1 , ALPH  , Y ! ) 
WRITE<6»6)PM»V.C11*ALPH, Vj 

WR I TE ( 6#  3 ) 

3 F0RMAT(2X,74HcXU,CXA»cZU#CZA.rZAD,CZTDj . CMU . CM Al » CMaD . CMTDl . CZTc  . r 

IMTC.CZBIiCMBI.CMI.CXI ) 

WR1TE(6#6>CXU,CXA,CZU,CZA,C7AP 
WRITE(6,6)CZTDl.CMU,CMAiiCMAD,CMTni 
WRITE (6i6)CZTC,CMTC.CZ01#CMB1,Ch! 

WR  I TE ( 6 » 6 ) CX ! 

CM  I »CMB1*CM I 
CZI*CZ81 
Y 1 1* Y l 
97  CONTINUE 

READ(5,1)FM,GV»C5Q.C5i.c52 

wRlTE(6.l)FM,GV,C50,Ci5i,C52 

fg«pm*c.v 

YIsYJi-FG*CMI«C50 
Cl2=FG*C52 
C13»FG*c51 
C14«0  , 

C15«FG*C50 

AsCXU 

daCZU 

CMTD»CMTD1*C13«CMI 
CZTD3CZTD1*C13*CZI 
C = CXA 

C WHEN  v=0  PMsM/Q* AD  QsRhO*vtId*vTTP  USE  CZW  ETC. 

IF(V-i. ) 9 » 9 » 0 

8 F=CZTD+PM 

GOTOlO 

9 FscZTD 
10  CONTINUE 

0=CZAd-PM 
ESC^ A*Cl4*CZ I 
CMA«CMAl*Ct  A«,CMI 
Al*pM*D*Y  1 /\l 

01*,,(A«D«Yl)'PM#DoCMTn/v*pK(,Yt#p/v*pM*F*nMAD/v 

ClsA«D«CMTD-Ao6itYl-3MoE#CMTn/V^PM#F*CMA/V-A»F#CMAri 

C1sCi*B*C*YNC12*PM*(CZI«CMaD-CmIcD)/V 

Dl*A«E»CMTU-A«F ttCMA  + CsF  « C M U 

D1sD1-B*>C«CMTD-B*CU1*CMaD 

DiaDl*D«CLl<>CMU 

Dl5Dl<»Cl2«P'1*{CMA*CZl-6*CMn/V^Cl2*A«(n*CMI-C'!AD*CZl  ) 

E 13« ( B«CMA«CL 1 ) *E *0  1«CmU 

EI^I  + CISoaMe^mI-cMa^cZTj  + CTZ*^*  CMU«CZ  i - B* CM  I ) 

C90  = CZA°CMU-CZY<*CMA 
C91=D*CMU-CZU*CMAD 
Elsfcl*C12«CXt<»C90 
D1  = U1*C12oCXI<*C95- 
Cl=Cl*C13»Cx!»C9l 


B-6 


D1*D1*C13*CX!«C90 

Al»Ai*C15*CZ!*CMAD*PM/V-Ci5<»Cy  !*D«PM/V 

81*«l  + Ci5*CZl<*(  (PM/V)<‘CMA-CXU«CMAn)+Cl5*CMl“<CXl;*D-CZA«PM/V)*-C15*r 

1XI*C91 

Cl*Cl*Clb*CZ!«(CMU*CXA-CMA«CXU)*Ci5»CMl»(CXU*CZA-CZU#CXA>*Cl5»CXI* 

1C90 

A2*l, 

82*dl/Al 
C2*Ci/Al 
D2*Dl/Al 
E2*E1/A1 
XCOF ( 1 ) *E2 
XCQF { 2 ) «D2 
XCOF ( 3 ) *C2 
XCOF ( 4 > 8B2 
XCOF { 5 ) ■ A 2 

5 F0«MaT(2X,43HQUaRTIC  cOcFf  j f;  J cNjS  ARE  jN  ORDER  A*B»C»D#E) 

WR  I T6  ( 6 • 5 ) 

WRITE(6i6)A,B,C,D.E 

WRlTE(6»6)Al,Bi.Cl»nl.El 

WRITE<6*6)A2,B2,C2,n2.E2 

6 F ORM*7 < !>El5 . 8 ) 

WRlTg<6#l)Ci2,Cl3.Ci4 

CAUL  ROUHT(xCOF,COF,H,RoOTR,RrOTI  . IER) 

WR 1 TE ( 6 » 7 ) 

7 F0RMAT(»X.9HREAL  ROOT , by  * 1QH ! M Afi , ROOT , 0X • 6HP6R I Oq. 9X , 12«T IMF  TO  u 
lALF) 

00  2 U ] al  # 4 

IF (ROOT  I ( 1 5 )8l, 80.81 
80  P3200000. 

GO  T 08  2 

01  P*2t«3.l4/R00TI<!> 

82  D* , 69/R00TR  < I ) 

20  WHITg<6#  U>ROOTR<  I > . «oOT  I C H.PiD 
READ<5»29)M0RE 

29  FORMAT(IS) 

1 F (MORE  5 98 » 30 • 97 

30  I ER*2 

11  F0RMAT(bX.6517.8) 

END 

SUBROUTINE  ROURT  ( XCOF  , COF  , M , Rr  OTR  . ROOT  J , j*-R> 

dimension  xcof  < io  ) » co? <io) , ro^tr ( io  > .root i t i o > 

I F I T a 0 
NSM 
I E R ® 0 

IF (XcOP  < N*1 ) j 1 0 » 25 1 1 0 
10  I F ( N ) l5 i 15 i 32 
15  IER*i 
20  RETURN 
25  I ER*4 
GOT020 
30  I E R 3 0 
GOT020 

32  l F ( N-36  >35.35,30 
35  NX=N 
NXXSN*1 
N2s  1 
KJ1=N*1 
D040U=1.KJ1 
MTsrlJl-U+1 
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40  COF(MT)aXCOF(L) 

45  XO*  . 00500101 
YQ*.  01000101 
IN»U 
50  XaXO 

XO*-10.*YO 

YQ»«10.*X 

X*XO 

Y*YO 
I Ns I N*1 
G0T059 
55  I F I T*1 
XPR«X 

ypr«y 

59  ICT*0 

60  UX*0 , 

(JY30 , 

V* 0 i 
YTsO, 

XT*l, 

0*COF(N*1) 

IF ( U ) 65> 130 » 65 
65  DO70l«l*N 
L*N-I*1 
XT2«X«XT-Y*YT 
YT2«X*YT*Y*XT 
UaU*COF(L)«XT2 
V» V + CqF  < U ) * YT2 
Fl*I 

UY*UY-FI*YT«COF(L) 

UXsUX*FI*XT*COF(l> 

XT  = XT2 
70  YTSYT2 

SUMSQsUX«UX*UY*UY 
IF  <5>UmSQ)75»110i75 
75  DXs( v«UY-U*JX>/SUMSn 
XsX*DX 

DYs-(u*UY*V«UX)/SUMSQ 

Ysy  + DY 

78  IF(ABS<DY>*ABS(Dx)-.OnOoOl jVon.80.80 
80  ICT»ICT-1 

IF(  I CT-i>00  ) 60  , * 85 

85  IFdFITUOO. *0.100 
90  l F ( l N"5 ) 50 » 95 . 95 
95  IE«»3 
G0T020 

100  D0105L«l.Nxx 

MTsKJi-Ln 

TEMPaXCOF(MT) 

xcof ( mt  > scor ( l ) 

105  COF ( L ) sTEMP 
I T6MPsN 

NsNX 

NX= I TEMP 

I F ( I F I T)  120 » 55*120 
110  IF(IF!T)u5f5o,ii5 
115  XsXPR 
VsYPH 

120  I F I T s0 

122  IF (ABS(v/X) -0.00001 >135, 125.155 


125  AIPHA«X*X 

SUMSQaX«X*Y*Y 
NiN-2 
QOT0140 
130  XaCi 
NX'-NX 
NXX«NXX«1 
135  Y3 0 

SUHSQaOi 

ALPHAaX 

NaN"l 

140  C0F<2>*CQF<2>*AlPHA#CoF(l) 

IF(N-2)147»i46i146 
147  NN«2 
G0T0145 
146  NN*N 
145  D01i>OL*2*NN 
M11»L*1 

150  CqF<m11>5CoF  « MH  > ♦ A.  PHA«C0r  (L  i -sUmSO^CoF  ( 1 -1) 
155  H00TJ(N2)*Y 
R00TH{N2)sX 
N2«N2*1 

IF(SUmSQ)160i165i16o 
160  Ys-Y 

SUMSQ«0i 

G0T0155 

l65  IF(N)2o#20'45 

return 

END 

LIST (STOP ) 

LIST 


« DATA 


,000194 

1, 

.0000079 

,00ft2 

-.0000042 

-.000009 

0| 

.000000078 

,000000008 

-.041 
, 0045 

,00000336 

0. 

0.  0. 

0. 

0.  0 

1 

1.  1- 
-1 

1. 

2.  4 

1.26 

100. 

• 405 

.000466 

, 05207 

-. 0001 73 

0. 

.00035 

• . 04?7 

-3,14 

.363 

.155 

.793 

A A 

0. 

•*.000043 

0. 

- . 0 0 1 3 


0. 

-.  793 
-.00008 
-.160 


0 


PROS 


’.  00051 
0. 

.00000052 

0. 


• 033 
0. 

-.13 

-.0144 
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appendix  c 

LONGITUDINAL  Tl  .TENT  MOTION  PROGRAM 


C-l 


OOOOO  O O O O OOOOO  OOO  OOO  O OrtO  O O OOOOOOO  ooo  ooooo  ooooooo  oo  • 


MAXO ( 2400 ) 

U I ST  < STAHT ) 

Long^tudInai  mOtIow 


VARlA8Lg  FWD  LOOP  G A I ivjS  LImIT  Gyu  <*VL  CORRESPOND  T0  «Tu  UPPER 
FEEDBACK  GAINS  C»0 < S2 > c5i(S)  C52 


1 

1-15 

PM 

1 

16-30 

V 

l 

31-45 

CL1 

1 

46-60 

ALPH 

1 

61-75 

Yl 

2 

1-15 

CXU 

2 

16-30 

CXA 

2 

31-45 

CZU 

2 

46*60 

CZA 

2 

61-75 

CZAD 

3 

1-15 

CZTD 

3 

16-30 

CMl 

3 

31-45 

CMA1 

3 

61-75 

cmtdi 

4 

1-15 

CZTC 

4 

16-30 

CMTC 

4 

31-45 

CZB1 

4 

46-6q 

CMBl 

4 

61-75 

CM  i 

4 A 

1-15 

CXJ 

5 

01-10 

RTU 

5 

16-30 

rtl 

3 

21-30 

GVU 

5 

31-40 

gvl 

6 

01-10 

fm 

6 

11-20 

C50 

0 

21-30 

C51 

6 

31-40 

C52 

6 

41-50 

SG 

mass  m»v/q*ad 

velocity  FT/SEC 

lift  coefficient  w/u*ad  ad  disc  arfa 

*NG».E  OF  ATTACK 

moment  of  inertia  i Yy/u*ad*Rotr 
stability  derivative  dcx/du 

dcx/dalpha 


FOR  STABILITY  DERIVATIVES  SEE  HBLin  4 


co'LfctivE  pitch  control  effect 

CO*  LECTI vE  Pitch 
cyclic  pitch 
CYCLIC  pitch 

HORIZONTAL  TAIL  r F^  ECT I VgNESS 
PITCH  CYCLIC  VARIABLE  GAIN 


AUTOPILOT  GAINS 


0 EnD,-1  nEw  PM,  1 NEw  LAIN 


rwC  VEL) 

xu*u 

ZU«U 

MU-U 


LONGjTUr>lNA|  EQUATIONS 

alpha  Theta 

X A j A -Cl.T  *MU°UD 

2A#A*ZaD*Ao  Ztd«Tf  «M*(WD-VoTD) 

MA*A*MAD<>A’’  Mtii»T~  * I YYeTD1' 

DETERMINANT  ;n  FORWARD  FI IG~T 


F WD  VEL, 
CXU-(PM/V)*S 
CZU 
CMU 


anGlS  of  attack 

CX  A 

CZA  ♦ ( C 2 A D " ^ o c 
CmA  *CMaDoS 


PITCH 

-CL 

(CZTO+pyloS 
CmTD*S  *Y[*5«S 


PMrM»V/ ( "*> Afn 
Yl  = lYY/(fj*ADt>p 
ALPHAsw/V 


C 

C 


Fh'D  VEL, 
CXU-PM«S 
CZU 

CMU 


DETErmInA.n*  For  mover  V = o 

vertica!  vEii  pitch 

CXU  CL 

czw  ♦(C^wn  pm > * r cztd«s 

CmW  ♦cMwDoS  CMTD*S-V i°s«s 


97  CONTINUE 

READ (5, 7 )Pm, V,CLl, A_Ph»  Yl 

READ(^,7)CXJ,CXa»CZj,cZa»C7aD 

READ;tj»7)CZTDl»CMU,CMil,CHAn*"MTDi 


PM=M/ ( QoAD ) 
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READ(5#7)CZTC,CMTC,C2B1,CMB1 # CM  J 
READ(5»7)CXI 

7 FORMAT ( 5E15 . 8 ) 

1 FORMAT ( 7F10 . 6 ) 

WR I TE ( 6» 100  ) 

100  F0RMAT(!>X#48HAlRCRAf-T  Qa  INPUTS  IN  ORDFR  ARE  PM , V , CL  1 » AlPH , Y J) 
WR1T6(6>7)PNiV»CL1»  Ai.PH.Tr 

WR1TE(6»101) 

101  FORMAT < 2X, 74HCXU, CX A » CZu.CZ A, rz AD, CZTDl.CMU.CMAl.CMAD.CMTDl.CZTC.r 

1mTCiCZB1iCmBi»CmI»CxI  ) 

WRITE(6.7)CXU,CXA,CZU,CZA,n7Ar 

MRITE<6.7)CZTDl.CMU,CMAi,CMAD.CMTnl 

WRITE{6.7)CZTC.CMTC,C78i,cmr1.CmS 

WRJTE<6»7)CX J 

CZTD»CZTD1 

CMA*CMA1 

cmto*cmtdi 

CMD6«CMBi*CMJ 

cZde«czbi 

TV*0 , 

V9*V 

98  REA0(5.1)RTJ,rTliGVj,bV!. 

READ(5.1)FM,C50»C51,C52,Sg 

WHITE(6»102)RTU.RTL,GvU,Gvl 

102  ForMAT(&X,4HRTUS|F9.3.1X.4H9T  ■ ,Fg.3» 1X.4MGVU* ,F9 . 3, 1X»4MGVI«  ,F9 . 3 
1) 

WRITE(6.103)FM,C50,C5l.C52 

103  FORMAT(»X,3HFM«,F9.4.lX,4wc:50s*F9.4,lX,4HC,Jl*.F9.4.lX,4HCtj2=.F9.4) 
V«y9 

001*0 i 

TDD2«0. 

GV»0. 

DD2*0, 

THO*0. 

THD*0, 

ALO*0, 

Vd*0, 

TM»0, 

DTM* , Q5 
CDaO , 

DL  aOTM/2 . 

DEMX«l,/57,3 

N*0 

TH*THO 
AO*TH 
VO  *V 
AL*AO 

C WHEN  V*0  PM sM /O i» A D OsRhO*vtIP«vTtP  USE  CZW  ETC. 

IF ( VO ) 8 » 9 1 8 

8 F 1 0 s V 0 / P n 

80T010 

9 FlO*l./RM 

10  continue 
fh=i>m-czad 

DO  80  1=1,800 
CDl'CI) 

!F(TM-,5)70, 7i, 72 

70  CD*TM/ , 5 
GOTO  75 

71  CD«1, 
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72  IF(TM-1, >73,74*74 

73  CD» ( 1 • "TM ) / . 5 
GOTO  75 

74  CD"Q i 

75  DE«CDl«*DEMX 

DA«AL-AO 

DV*V-VO 

oth*th-tho 

RT*ABS<  ( DD2"DD1 ) ) 

IF(HT-RTU)201#200.200 

200  GV»GVU 
GnTy  204 

201  I F < KT.rTU  >202, 202,203 

202  G V*GVL 
GOTO  204 

203  UV*1»VL*<  (RT-RTL)/(RTU-RTL)  >«<GV u'GVl.  > 

204  CONTINUE 

IF  { TM-1 , )35q.  35.1,35i 

350  GV=0. 

351  CONTINUE 

fg»fm*gv 

VDl*<CXU«|)V  + CXA*DA“Ctl*OTH*f‘G<*CX!*<Ci>0*Tnn2  + C5l*TMD- 
ALDl*{CzU*Dv*CzA»DA*.(CZTD  + PH*r5l*rG*CznF)  (lTHD  + C52#F( 
1*(DE*c50*Tdd2*FQ) )/Fli 

TDDl«(CMU*DV,CMA»DA*CHAD*Al-nl*(CMTi34C5l  *FG«CMDE ) *tH[ 

1FG«TH)*TV)/(YI«C50»fO#CmDf> 

DD1«Fm*GV*<C50*TDD2*C51«Tmd*CS2*Th) 

DD15s aBS ( DD1 ) 

I F C DD15- • 26)41,40, 4q 

40  DDla,26*DDl/Dni5 

voiB  (Cxu*Dv+CxA»DA-r,ii*nTH*nx!*(Dnl*DE  n*?  io 

ALDls(CZU«DV*CZA*DA4.CZ1D‘:Tnn*f:ZDE«<0Dl*DF  ))/Fn 
TDDla(CMU«*0v*CMA*0A*CMT0*TwO*('MDE»(ODi*DE  ) )/Y  I 

41  CONTINUE 
Vi»V*VDl<*DL 
ALla AL*ALD1*DI 
TDlaTHD*TDDl*Dl 
Tj>«TH*THD*DL*,5*TDDl*0L«Dl 
V2*V14VD1«0L 

AL23AL1* AL0l*DL 
TD2*TD1*TDD1*DL 
T2*Tl*THD*DL* , 5*TDDl*nU«D(. 

DAaAi.2-AO 
DV» V2-V0 
DTHaT2-TH0 
DEaCD«DEHX 

DDlaFM*UVMC50*TDDl<.C51«Tn>i.  42*Tt ) 

RTa ABS( (DD1-DD2)/DTM) 

!F(HT-RTU>301,300.300 

300  GVsGVU 
GOTO  304 

301  IF<RT-RTl.)302,302,303 

302  GV  = GVt. 

GOTO  304 

303  GV=GVL*(  (RT-RTL)/(RTU-RTU)  )*(GVU-GVi.  > 

304  CONTINUE 

I F ( Tm- 1 j >2.2.6 
2 GV=0, 

6 CON T I muE 
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GOTO  75 

72  IF(TM-1, >73,74.74 

73  CD* ( 1 « "TM ) / . 5 
GOTO  75 

74  CD"0  • 

75  DE*CDl*DEMX 
DA«AL-AO 
DV*V-VO 
OTH*TH-THO 

«T*ABs<<DD2-DD1)/DTm) 

IF<HT-RTU)201#200.200 

200  GV«GVU 
GnTo  204 

201  IF(HT-RTI.)202, 202,203 

202  GV«GVL 
GOTO  204 

203  GV*Gvi*<  (RT-RTL)/(RTU-RTL))*(RVu-GVi.  i 

204  CONTINUE 

IF(TM-1.  )35o.35.i,35i 

350  GV*0. 

351  CONTINUE 
fg«fh*gv 

VDl*(CXU»|)V*CXA*DA"CLl*r)TH*rG*CXl«<c30*Tnr)2  + C5l*TwD  + C52*TWn«Fin 

ALDj_a  (C2U*Dv*CzA«DA*(  CZTD  + P«*r51(»FG#C2PF  )„THD  + C52#FG»CzPE<>Tu  ♦CzHc 

1«<DE*c50*TDD2*FG)  )/Fll 

TDDl«(CMU*DV,CMA*DA*CMAO*At- nl*( CmtD+C5i*FG«CMDE ) *tHD*CMDE*(DF,C52* 
1FG*TN)*Tv)/(  YI-C50»fO*CmDf) 

DD1»Fm*GV*(C50*TDD2*C51«TmD4.C52#Tm) 

DD15* ABS ( DDl ) 

I F < DD15- . 26)4l, 40,4q 

40  DDx*.26*DD1/DD15 

VDls(CXU#Dv  + CxA*DA«rA'l*PTM4'nXt*(Dnl*DE))*';lO 

ALDls(C2U*DV*CZA«DA*CZTD*T«n*CZDE»<DDl*DF )|/Fn 
TDDXa(CMu*Dv*CMA«OA*CMTn»Twn+nMnE«(ODl*DE ) ) / Y I 

41  CONTINUE 
Vi*V*VDl*Dl 
ALla AL*ALD1*DL 
TD1*THD*TD01«DL 
Ti«TH*THD*DL* . 5#TDDl*0U«D| 

V2*V'i*VDl«0L 

AL28AL1*ALD1*DL 

TD2*T01*TDDl*DL 

T2«T1*THD*DL*.5#TDDi*0L*D|. 

DAaAl.2"A0 

DVaV2-V0 

DTHaT2-TH0 

DE3CD*DEMX 

DDl3FM*GV«(C50*TDDl*C51«TDi . S2«Tt ) 

HT*ARS< (D01-DD2)/DTM) 

IF(RT-RTU>301»30S.300 

300  GVsGVU 

goto  304 

301  IF(HT-RTL)302,302,3o3 

302  U V s G V l. 

GOTO  304 

30  3 GV=GVL  + ( (RT-RTL)/(RTU-RTU)  )*(GVll-GVL  ) 

304  CONTINUE 

IF (TM-1j ) 2 » 2 i 6 
2 GV=0, 

6 CONT  I 'JUb 


FG«FM#GV 

VD2»<CXU*DV*CXA#DA-cLi*nTH*FG*Cxi*<C52»DTH*-C5l*TH0+C5n*TD0l) >*Fi  0 
ALD2«<CZU*DV*CZA*DA*(cZTD+PM*r5i*FG*CZnE)<*TL)2*C52#FG«CZ0E«T2*CZnF« 
l<DE*Cf50*TDDl*FG))/Fil 

TDD2»<CMU#Dv*CMA«DA*CMAn»Ain2*CCMTD+C5l «Fo*CMDE ) «TD2*CMDE* ( DE*C52* 
lFG*T2)*TV)/( VI-CSO*rG»CMDF) 

DD2aFM*GV#<C50*TDDi*C5l*TMn*C52*DTH) 

DDlt»»ABS  < DD2 ) 

IF<DDl5-.26)43#A2#42 

42  DD2«,26»DD2/DD15 

VD2»(CXU«DV*CXA*DA-cLl«DTH*rXT*(DD2*DE) >«FlO 
ALD2* (CZU*DV*CZA»DA*CzTD*Trt2+rZnE»(DD2+DF ) ) / f 1 1 
TDD2*(CMU«OV*CMA«DA+CMAO*ALn2+CMTD*TD2*CMUb<»(OD2*DE) )/YI 

43  CONTINUE 
V*V1*VD2*DL 
Al*AL1*aLD2*Dl 
THD*TD1*TDD2*DL 
TH»Tl>TDi«DL*,5*TDD2**nL«DL 
THP*57,3*TH 

l F ( VO ) H » 12  # 11 

11  AlP*57  # «J*Al» 

gp«thp-aup 

G0T013 

12  ALP»0. 
gp«tmp-alp 
alp»al 

13  continue 
GLaA3S<UP> 
lF<Gl.-90.  >79,81.81 

79  CONTINUE 

IF  ( fy-1  >20,20,21 

20  N>N«1 
GOTO  80 

21  N*0 

DD2*FM#Gv*<C50«TDD2*C5l*Twn*C52#Tw> 

DXa!>7 , 3#DD2 

0X1«57,3*0D1 

I F ( V >606.607 , 606 

6 06  WRITE(6,104>TM,V#DX1,nX,A|.p,TuP,GP 
GOT0608 

60 7 WRITS(6.60b>TM#V,DXl,OX, AlP.ThP.GP 

608  CONTINUE 

104  F0RMAT(3X,3hTm*,F6,2,iX.2hV5,F5.1.1X»4pd*‘'s'F9.5,iX,4h  dX=,F9.5, 1 x 
1#4nAI.P8,F10i6,1x,4htHPs  ,Fm  6 , lX,3H6Ps  .Ft!)  .8  > 

605  F0HMAT(3X,3HTMs,F6,2#lx,2HVs.r5.1.1X,4H0Xl=,F9.5,lX,4H  DX=,F9.5.«  X 
1 , 4H  ws ,F1Q  .6,1X,4hTHp=.F1o , 6 , IX . 3HTH*  . Flfl , o > 

80  Tm»Tm*DTm 

81  IF ( SG ) 97 » 30  » 98 
30  CONTINUE 

END 

• (.  1ST  < STOP } 

LIST 

• DATA 


,000194 

0, 

.0062 

0. 

.00051 

-.0000042 

,0000079 

’.000009 

-.000043 

0- 

0, 

,000000078 

,000000008 

0. 

-.0000043 

-.041 

,00000336 

0. 

-.0013 

o. 

, 0045 

10,  1.  0.  0. 


C-5 
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APPENDIX  D 

PROGRAM  FOR  SOLUTION  OF 
LATERAL  EQUATIONS 


D-l 


4 SA5  CONTRIBUTIONS  4 SAS  CONTRIBUTIONS  I BASE  VALUES 


_ m' 

w * 


Ct 

•vr*,a 

c*?4* 


(-m*Cy^  |» 

<V 

cV^»'rZ*2 


■frxr2  Cyv  i'x i'i 


■yU'z  cl4+Ixc,?4l 


Cgi 

(D^ 

C*j  (I)**'" 

r (21**- 

C 7 

c* 

.(51**- 

(S)*'"’ 

(6)-- 

1 r 

■CYv*X  c^8tr 

TSC”®tr 

•TC^CH 

C>?v1xcYeIfi 


■Cyu|I'zC|^>  lx  C,^)  Cy^C^C^-Cj^C^) 

CluI'zCY4  C»jy(Cy^Ct^-Cy^C|^) 

Cijw  I'x  Cy^ 

-mCijuI'x 


Ci^Cl  Ct  Cl^ 

-CjU*cL 


-CywI2  C^j 

7c^cUi 

"v"  Ci'#'C^AI 

<*«  I‘z  CyA, 


!6) 

(8)< 

C’«'  Cl  Ci9ir 


cYuct4 

*cyvcij^,  cieTR 
Cl„c^  CyeiR 

*C*«Cy4  Cl?eTR 
'C,7uCl4  cYgIR 


*NOTE  TRIANGULAR  RELATIONSHIP 
9 9 9 

ID  , , , 

(!)  CVINTIC 


m'Ci;uCiAl 

CYvC^C1a| 

Cy«Ct^,C>jA, 
Clu  Cy^,Cj}a1 
-m'Ctw  Ci>AI 
Qlu  ClJ^rCyA1 
C^Cj^CyA| 

C^Yy'-tAI 


T.'BiA  D.l  expansion  of  Laferoi  - Directional  Derivatives 


n-2 


• 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LIST (START ) 


5?»  39? 

SOLUTION  Of  LAtFRA'i'  fQijATIONS 


SV 

-PM(S)*CVV 

CLV 

CNV 


SY 

(-pm*v+CySd> (S; 
Cl  Sq  ( S ) 
-Zl  <S**2)*OSD<S> 


CA«D 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

4 
4 

4 

5 
5 
5 
5 
5 


column 

01-15 

11-20 

21-30 

31-40 

41-50 

1-10 

11-20 

21-30 

31-40 

41-50 

5i"6q 

61-70 

1-10 

11-20 

21*30 

31-40 

41-50 

1-10 

11-20 

21-30 

1-10 

11-20 

21-30 

31-40 

41-50 


SYMBOL 

PM 

V 

CL1 

XI 

ZI 

CYV 

CLV 

CNV 

CYSD 

CLSD 

CNSD 

CYPHD 

CLPHD 

cnphd 

CYDS 
Ct.DS 
CHDs 
CIDr 
CL  DR 
CN  OR 
FM 

ov 

C6C 

C61 

C62 


VELOCITY 

LIFT  COFFFlClFNT 

moment  of  Inertia 
MOMENT  of  inertia 


C YTh ( S ) +CL 
.XI (S*#2)+ClTd(S> 
CNTd^S) 


FT/SEC 

W/Q*AD  AD  DISC  AREA 

!XX/Q*AD*ROTR 

I22/Q*AD*RQTR 


MODEi  GAIN 

Gain 

autopilot  coefficient  S2 


c 

typical  inputs 

Cl. 26 

100. 

.405 

.0071 

.0278 

C-. 00094 

-.000072 

. 0004 

-'.Oil 

0. 

co. 

C.114 

-.000085 

,025 

.000015 

-.114 

-’,4  05 

-.116 

co. 

Cl 

0. 

0, 

O'. 

o. 

C*, 

C-l 

1. 

1. 

2. 

4, 

C. 000194 

1. 

.0062 

.00011 

.00043 

C-. 0000^1 

-.00000076.0000052 

.000148 

0. 

CO. 

f*  n n 4 ii 
V » V W A V 

.0000018 
. 0004 

. 0000003l 
-.0018 

- 0 0 45 

-',0013 

Co. 

cl 

o. 

o. 

o’. 

o, 

C«. 

co 

1. 

1. 

2'. 

4, 

193 

.0201 


-.0001A8 

-.00023 


PROB 

DIMENSION  XCOF(ll) ,cUF(il)  ,rocTR<1 0),ROOTI (10) ,Z<4) 

COMMON  XCOF.COF,M,RoOtR,RoOTI , ier 

11  FORMAT ( 2X i o5l7  t 8 ) 

12  FQRMA  T ( 1 X # 5QH  REA.  IMAG  PERIOD  T 1/2 

11 0 H DR  /IX . 6Fi 0 . 4 > 


WN 


o o o o 


97  CONTINUE 
J*  0 

READ (5, 1 > PM, V#  CL*  XI , Zj 
READ(5*l)CYV*CLV*CNv*CYSD,n.Sn,CNsn 
READ(5,l)CYoHD,CLPHD«CNpHD,rYrS,CLDS,CNiDS 
READ(5,1)CYDR,CLDR,cNdR 
WRlTE<6»95)pM,V.CL,xJ,Zl 
WRlTE(6,94)CYv,ClV,CNv,CYsn,Ci.  SD.CNSD 
wRITg(6*93)CVPHD*CLaHD»CNPHn»CYDS.CLDS.CMDS 
WHITE(6,92)CYDR,CLDR,rNDR 

92  F0RMATaX,5Hr,YDRs»El5.8,lX,'?HrLDR=*El5.8,  i X . 5wCNDR=  • El  5 . 8 ) 

93  F0RMAT(lX»6HCYPHDs,El5.8*lX.6MCLPHD=*El5.«.lX,6HCNPHD=.Ei5.8,lX,^u 
lCYDS«,El5,8,lX,5HCLnS  = »El5‘.8,lX,5MCNDS=,El5.8) 

94  F0Rmat<1X,4hCYV=*E15.8*iX,4hCi  V«  , e!5 8 , IX  , 4hCNV=  » gl5 , 8 , lX  , 5hcySd=  • 

1615,8  * IX, 5HCLSD, E15 , 8 , IX  * 5HCNSD«  »Pl5  * 8 ) 

95  F0RMAT(lX,3HPM=,El5.8,ly,2MVB,Ei5.8,lX,3HCL:,El5.fi,lX,3HXI=,El5.fl, 
HX , 3HZI 3 * E15  • 8 ) 

IF  INPUTS  CYB  IF  NOT  RgMovc  NEXT  3 CARDS  IE  CYV=CYV/V 
CYV=CY V/V 
CLV»CLV/V 
CNV«CNV/V 

98  READ(5,l)FM,OV,C60,C6i*c62 
WRITE(6*91)FM,Gv#C60#C61*C6? 

9i  FORMAT(lX,3HFMs,F6,3»lX,3HGvs>r6.3»lX»4HC60=*F6.3,lX*4HC6ie,F6.3,l 
IX i 4HC62* • F6 , 3 ) 

FG*FM*GV 
M»4 

ClOSaPG*C6Q 
ClOaFG*C61 
Cll«FG*C62 

C FOR  HOVfcR  S=T  V=1  and  PMaM/(0«An>  0«Rh0*VT1P*VT1p 
IF  tV-l, >3*3.4 

3 FlOsO, 

Fll*0 , 

F12*  0 , 

WR I TE  < 6 , 6 ) 

6 F0RMAT(5HH0VER) 

G0T05 

4 FiqsCNV*XI*CYSD 
FnaCNV*CLPHD«PM»CL  V»PMoCnPwD 
F12*PM*CNV,CL0S-PM*CLV*CNDS 

5 CONTINUE 
As«pm*xNZI/V 

B = CYV*XI*ZI +PM«(ClPhD»Z! ♦cNSD«X!  )/V 

C*-CYV* (CLPHD»ZI *XI oCNSd)-Pm»( CLPtJD«CNSD-CNPHD*CLSD)/v+CLV«CVPHn»7 
1 1 ♦F 10 

CsC-CNV*Xl®aM 

DxCYV*ICLPHD*cMSD-CNP^D«CLSn)+Cl.  Vo(Zl*CL»r,  YPHO«CMSD  + CNPuD*CVSD)»Cf' 
lV*(CYpHD»CLsD-CLpMD*CYSD) 

D=D*r 11 

Es-ClV*CL*CNSd+CNV«CL»ClSd 

XlsCVV*(CLDS*CNSD-CNDS0CtSD) ♦rLV°< “C,YD?*CNSU*CYsDoCnDs  J+CyVM'IYDs® 
1CLSD-CySD*ClDS) 

XlsXl+Fl2 

X2  = *cyY<‘CLQS»ZI-Pm®<C;  DS#cNSD-Clsn*cNDR  )/V  + CLV«CYDs<,z  I 

X35Pi-i«CLDS»7  J / V 
Y I sPM»x i »CNDR/V 

Y2s--YV(>Xl«QNoR-<pM/V)#(clpMr)0cNDR"Cl-DRttCN'PHD)'fCNv<>CYDPo)<^ 

Y3»CY7* (CLPHDoCNDR-cLnR»CNPwn)*CLv«(CYPR«rNPHD-CYPHD®rNDR  )+CNv« 
1PHD«C'.Dr-CYDr»ClPHD) 


D-4 


Y4»CNV*CL*ClDR-CIV*CUCnDr 

E1*CU*X1 

D1*C11*X2+C10#X1 

C1»CU*X3*C10*X2*C1qS«Xi 

81«C10*X3*C10S*X2 

Al*ClOS*X3 

A»A*Al 

C*C*C1 
D»D*Dt 
E*E*E  t 
A1®1 , 

B1«B/A 
Cl«C/A 
UlaD/A 
El»6/A 
XCOF ( 1 ) *E1 
XCOF(2)*Dl 
XCOF ( 3 ) *Cl 
XCOF { 4 ) »B1 
XCOF { 5 ) * Al 

HKlTg(6il>U,Bl.Cl»nl.El 

CAI.lR0l.HT  ( XCOF i COF  » m*  POOTR, ROOT ! i tER) 

00201*1*4 

:f<HOOtI ( I ) )8l»80*8t 

80  Pa203000, 

G0T082 

81  P*2.*3.14l7/R00TI  ( 1 ) 

82  D* , &9/R00TR ( J ) 

WN« (ROOT  I ( I )*«2*R00TR( 

DR*-ROOTR< I )/UN 

20  wRITE(6*v>0)ROOTR(  l > . RoOt I { J ) * » * 0*  wN . DR 

9 o FORMAT(lX,6MROOTRs#ci2,5*1x,6KRnOTls»El2.5.lX.7HPERlOn«<El2.5.lX.t 
14h  TIME  TO  HALF=iE12.5*iX,3hWN«,Ei2.5,iX,XhdRs*E12.5> 

1 F ORMAT ( 7F10 . 6 ) 

2 FORMAT ( 1 10 ) 

RE AD<  5 * 2 ) MORE 

I F ( MORE  >97*30,98 
30  CONTINUE 
END 

subroutine  RoiRTjxcoF.coF.M.RrCLTR.RoOTi ,ier> 

DIMENSION  XCOF(lO),cOF(lO),ROCTR(lO)iRnOTI (10) 

l F I Ta 0 

NsM 

IER*0 

IF(XCOF<N  + l)  )1Q. 25,10 
10  IF (N) 15* 15*  32 
15  IER»1 
20  Return 
25  IER34 
GOT020 
30  lERs0 
GOTO20 

32  IF(n-36)35* 35,30 
35  NXaN 
NXXaN+1 
N2*l 
KJ1=N*1 
D040i.  = l,Kjl 
MTskJ1-l»1 


D-5 


40  COF l MT ) "XCOF ( L ) 

45  X0«. 00500101 
YO», 01000101 
IN«0 
50  X*XO 

XO«-iO,«YO 

YO«-10.#X 

XaXO 

Y«YO 

I N* I N*1 

QOT 

55  JFIYbI 
XPR«X 
YPR*  Y 

59  ICT*0 

60  UX«0. 

UY»0, 

V*0. 

YT*0, 

XT«1, 

U«COF(N*l) 

I F <0 ) 65 * 130  » *5 
65  0070!«i»N 
L«N-J*i 

XT2*X*XT-Y*Y-f 

YT2«X*YT*Y«XT 

U*U*C0F(U*XT2 

V*W*cOf<L}<*yT2 

FI*! 

uy«uy.fi«yt.cof (l) 

UXaOX>FI*XT*CcF(|.) 

XTaXT2 
70  YTa YT2 

SUMSQaUX«UX*UY«UY 

tF(SUMSQ)75iH0.75 
75  DXa{  v*UY-0<*JX)/SUMS'3 
XaX+DX 

DYs-(U*UY*V«UX)/SUMSQ 

Y*Y*DY 

78  lF(ABS(DY)*A8S(0x)"-5f>000j  p 0Oi80.80 
80  ICTUCT-l 

IF(lCT-300>60» 85*85 
85  I F ( IF  I T ) 100  # 9Q  * 100 
90  IF ( I N*5 ) 50 , 95 , 95 
95  IEr=3 
GOTO20 

100  D0105l=1,NXX 
MTaKJi-Ln 
7EmP*xC0F (MT) 

XCOF(HT)scOF(L) 

105  COF (L ) aTfcMP 
JTEMPsN 
N = NX 

NX  = i TEMP 

I F ( I F l r ) 120  » 55,120 
110  IFCF!T)ii5.50.1l5 
115  X=XP9 
Y = YPR 

120  I F I T s 0 

122  JF(ABs(Y/X>-0, 0 0001  )1.' ?5,  t?5,l?5 


125  ALPHA*X*X 

SUMSQ*X»X+Y*Y 
N«N-2 
QQTQ140 
130  X«Q 
NXsNX 
NXX«NXX-1 
135  Y«Q 

SUMSQaO, 

ALPHAaX 

N*N»1 

140  C0F<2)«C0F<2)*ALPHA*C0F(l) 

IF(N-2)147»146»146 
147  NN»2 

QOT0145 
146  NNaN 
145  D0150L*2»NN 
H11*L*1 

150  COF(Mll>sCOFtMll>*A_PMA*Cora>-SUMSO*CrF<i.-l) 
155  ROOT  I < N2  ) a Y 
ROOTR { N2 ) aX 
N2«N2+1 

IF(SUMSQ>160,165,160 
160  Ya-Y 

SUHS(3*0  • 

00T0155 

165  IF(N)20,20,*5 
RETURN 
END 

• LIST(STOP) 

LIST 

• DATA 


1,26 

100: 

.405 

• 0071 

, 0278 

-,00094 

-.000072 

. 0004 

•.nil 

0. 

-.193 

0. 

» , OOUO  b5 

,000015 

• , 4(l5 

-.116 

- .02nl 

.114 

0. 

4 

,025 

0. 

-.114 

0. 

0. 

0. 

X 

4, 

1. 

1. 

2. 

4. 

-l 

,000194 

1. 

. 0062 

*00011 

.00043 

-.000061 

00000076 

. ooooo5*^ 

•000148 

o. 

-.000148 

0, 

.0000018 

.00000031 

•.0045 

-.0033 

-.00023 

.0018 

0. 

4 

. 0004 
0. 

-.0018 

0. 

0. 

0. 

X 

4. 

0 

1* 

1. 

2. 

4 . 

PROS 
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APPENDIX  E 

LATERAL  TP.ANSIENT  MOTION  PROGRAM 


• 

s 

MAXO(2400 ) 

. - . 

• 

LIST (START ) 

c 

LATERaL-DiRpctICNaL  motion 

c 

FcEoBaCk  GAINS  r60  C61  ”62 

C CARD  column  symbol 

c 

1 

01-15 

PM 

c 

1 

11-20 

v vfloCity  ft/s&c 

c 

1 

21-30 

CLl  LIFT  COFFP 1 C I FNT  W/Q«AD  An  RISC  ARFA 

c 

1 

31-40 

y i momfnt  of  Inertia  ixx/o*ad*rotr 

c 

1 

41-50 

/i  moment  op  Inertia  I2//q»ad*rotr 

c 

2 

1-10 

cyy 

c 

2 

11-20 

CLV 

c 

2 

21-30 

CNV 

c 

2 

31-40 

CYSD 

c 

2 

41-50 

CLSD 

c 

2 

51-60 

CNSO 

c 

2 

61-70 

CyPHD 

c 

3 

1-10 

CLPHD 

c 

3 

11-20 

CNPHD 

c 

3 

21-30 

CYDS 

c 

3 

31-40 

CLDS 

c 

3 

41-50 

CNDS 

c 

4 

1-10 

CYDR 

c 

4 

11-20 

CLOR 

c 

4 

21-30 

CNDR 

c 

5 

01-15 

RTU 

c 

5 

16-30 

«tl 

c 

5 

31-45 

GVU 

c 

5 

46-60 

GVL 

c 

01-11 

C60 

c 

6 

16-31 

C61 

c 

6 

31-45 

C62 

c 

6 

46-60 

SG  o END..1  new  PH,  1 NFw  GAIN 

c 

c 

typical  Inputs 

u 

i Cl. 26 

lOO, 

.405  .OOTi  .0278 

C-, 00094 

-.000072  .0004  -.Oil  0,  -.193 

1 co. 

-.  000085  . 000015  -‘.405  -.11  6 - • 0?0i 

C.114 

.025 

-.114 

CIO. 

.1 

4,  .5 

C,25 

1. 

2.  4.  1, 

CIO. 

.1 

4,  .5 

Cl.O 

1. 

2,  4’.  -1. 

c. 000194 

0. 

.006?  .onou  .00043 

C* .000061 

-.00000076.0000052  .OftO’48  0.  --00Q14B 

co. 

.0000018  . 00000q3i  -’.  0045  -.  0013  *-00023 

c.oo 

1 0 

. 0004 

-,00l8 

1 r-  n n 

1 

4 , ,p 

j C . 25 

1. 

2.  4.  1. 

CIO. 

.1 

4.  .5 

| Cl . 

1. 

2,  4,  0, 

' 

! c 

PtfOd 

! 100  CO\T!wUt 

R£AD<5,1)PM 

.V.CL.XI.Zl 

H E A U ( p , DCYVl.CLVl.^Nvl.CYsi.cLsD.CNSD 

I 

HE  A 

J .;5, 1 ICvTD,  CLTD,  CNtD.  Cyns  - rLns.  CnDS 

E-2 

f 

4? 

O O O t; 


READ(5#l)CYDR(CLDR,CNnR 

1 FORMAT (7F10 , 6 ) 
fRlTE<6#2)PM|V»CL*Mi7l 
WRlTE$6»3)CVVl,CLVl,CNVl#CYSD.ClSn*CNSn 
WRITE<6,4)CYTD,ClTD,CNTn»nvnS,ClDS,CNDS 
WRITE<6,5)CYDR,CIDR,CnDr 

2 FORMAT* lX.3HPM«,F6,3»lX,2HV=.':-5,l,lX.3HCLs»fr6.3.1X*3Hx!*,r7.4.lx.T 
1HZI*,F7,4) 

3 F0rMAT{1X,5HCYV1«,Fi0 .6,iX,5HrLVl=iFl0.6,iX.5HCNVis,Fio.6,lX,5HCVs 

1D**F10.6»1x»5hCLSD*.Fi0.6,1x»5HcNSD=iF10.6) 

4 FORMATUXi5HCYTD"*FiO .6,lx,5HCLTD=»Fi0.6. lX.5HCNTD*»Fio.6.1X,5HCYn 
iS^FlQ.ft.iX.SHCLDSo.FiO^.ix.sHcNnSs.FjO.b) 

5 FORMAT(lX,5HCYDR«.Fl0.6,ix,SHrLORs*FiO,6, tX,5HCNDR«,F10.6) 

9b  REAO<!jil>RTj,RTL.GVj#fiVL 

READ(5»l>FM.C6t,.C6l,C62.S5 

WR!TE<6.9)RTUiRTL.QVU,GVt 

WRITE(6i6)FM>C60»C6l»r62 

9  F0«MAT(lX,4HRTUs,F6.3.1X»4HRT.  • , Fft , 3# IX , 4hG VUs , Fft . 3 . IX , 4«GVl» . F6 . 3 

li 

6 FORMAT  ax ,4H  FMsfF6.3,lX,4HC60«,F6,3ilX,4-C6ls.F6.3,iX,4HC62«.F6.3 
1) 

WR I Tg (6 , 87 ) 

8 1 FORMAT(3X,2HTM,lOX,?HsV.lcX,2wSY»tOXi2NTH. 10X,2HDR*9X.3HDXR.9X.3HD 

1XA 1 9X 1 4HALAT  ) 

lc  :NPuT  CYv  REMOVE  NPXT  3 rAPDS  OKE  FOR  CYB 

CNV»CNVl/V 

CLV«CLVi/V 

CYV-CYV1/V 

CYV*CYV1 

CNV«CNV1 

CLV«CLV1 

SVD2»0. 

DD1"0. 

DD2»0. 
b 1 2« 0 . 

GRmX » 1 • /&7 , 3 
TDD2.0, 

SV*0, 

viY*0, 

TH«0, 

TDD*0. 

SVD«0. 

SD*  0 • 

TD*0, 

\'N*0 

TM«0, 

DTm*,05 
CD*U , 

DLsDTM/2, 

Drb0, 

FOR  HOytR  WHEN  ys0  Pm=M/<  0<>A<'  ) 0*RHO«»vT  I Po  VT  lp 

:rmio, iQ,u 

10  FJO*C. 

FU»i,/PM 

G0TO12 

11  FlO*PM 

y 11s  V /Pm 

12  CONT I Nut 
DObO  J J - 1 , dOO 
CD1*CU 


C 


T*TM 

XF  < T- . 5 ) 70 . 71 # 72 

70  CD«T/.5 
G0T075 

71  CO«l, 

Q0T07S 

72  IF(T-1. >73»74.74 

73  CD*  (j,  i “T ) / , 5 
Q0T075 

74  CD»0, 

75  0R»C01*DRMX 
RTsABS< <DD2-DDl)/DTy> 

IF(HT-RTU)201. 201.200 

200  GV*GVU 
UOT0204 

201  IF(rtT-RTl)202.202#203 

202  QV3liVL 
GOT0204 

203  GVsGVL*(  <RT-«TL)/(RTU-RTl.)  U(GVlJ-GVL> 

204  CONTINUE 

IF  CTm-1* )350*350»35i 

350  GV«0, 

351  CONTINUE 

fg»fm»gv 

SDDi<CNV*Sv>CNSD«SD*<CNtD*PR*C6*  »CNDS)*TD*CNDS*FG»C62«TH*CNns#rG*r 
160«TDD*CNDR«Cr)/2I 

SVDla(CL*SlN<TH)*CYvtt<?V*CY5n*GD*(rYTD4PG»'61*CYDS)*TD*C60*FG*CVOG» 

lTDD*C62*PG*CYDS*TH*''yitR(>DR.rift#sO)»rtl 

TDD«<CLV«SV*CLSD«SD.(rLTD4rr,*C61*ClDS)*TD*PG«C62*CLDS*TH«.CLDR.DR)/ 

l(XI-FG*C60ttCUDS) 

DD1*PM*GV*(C60*TDD*c61*tD*C62#TH) 

DDl3»ABS<DDl) 

IF(DDl5".2ft)4l, 40.40 

40  Pl)l».26»DDl/DDl5 

SDD«(CNV*SV*CNSD*SD*CNTn*Tn*CNDS*nOl*CNDR«DR)/Z! 

SVDla<  CL*SI  N(  TH)*Cyv*SV*Cvsn«SD*CvTO*TlHC^D5#DDl*rYDR*DR-,r10*SD?  )<• 
1F11 

TDQa(CLV*5V*C|_SD*SD*C1  Tn«TD*Ct  Dc,*nDl*Ci  DR«I)H)/Xl 

41  CONTINUE 
SDl»SD*S,'D*Dt 
TDlaTD*TUD*DU 
SvlaSv*Sv01*DL+DI2 
SlaSY*S01*l)L 
TiaTH*TDl«*Dl 

TD2»TD1*T0D*0l 

Sv2aSv*SVDl<*DL*D12 

T2sTi*TD1<»Dl 

SD2»SD1*SDD*DL 

RT  = ABS(  (DDi-DD2)/OTvi) 

IP  <RT-RTU)3Q1  , 500,300 

300  GV=GVU 
GOTO3O4 

301  IF(HT-RTL)302.302,303 

302  GysGVl 
GOT0304 

303  GV*GVL*( i«T-RTL)/(RTU.RTL) >»(GVu-GVl  ) 

304  CONTINUE 

IFITM-l, 5352.352,353 

352  GVsQ, 

353  continue 

E-4 


k 


DR«CD*DHMX 

fgbfm*gv 

SDD2»(CNV«SV2+CNSD*sD2*jCNTn*FG«C6l*CNDS)«TD2fCNDs*rG*C62»T2*rNDs» 
1FG*C60*TDD*CNDR*DR) /Z! 

SvD2a(CL*SlN{T2)*CYv«SV2*CYSD«SD2^<CYTr*FG*C61*CYDS)»TD2*C60*FQ*CY 
1DS*TDD*C62*FG*CYDS*TH  ♦CYnR»DR-F10«SD2)»Fll 
TDD2«<CLV*Sv2*CLSD*SD2*(CLrn*FG»C61*CLDS)*T02*FG*C62«CLDS*T2*CLDR 
l*DR)/(Xl-FG#C60*CtDS) 

OD2«FM»UV»(C60*TDD2+C61*Td2*C62*T2> 
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TYPICAL  LATERAL  STABILITY  DERIVATIVES 
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TABLE  F.l  (cont'd) 

TYPICAL  LONGITUDINAL  STABILITY  DERIVATIVES 
AND  NOTATION  USED  IN  EQUATIONS  OF  MOTION 
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CZTD  Var.  of  Lift  with  Pitch  Velocity 
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EXAMPLE  III  - SUMMARY  OF  STATIC  STABILITY (cont’d) 
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100  169  /5.7  7.4  -.031  -64.3  3.0  1.98  -.093 
125  212  3. 7/8.0  8.4  -.31  -11.4  .46  3.54  -.143 
150  254  4.4/10.9  9.2  -.61  -9.2  .34  5.6  -.207 


TABLE  F.5 

EXAMPLE  III  - STATIC  STABILITY  (Continued) 


STHADY  TURNS 
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NO  AUTOPILOT 


EXAMPLE  1 

v-:oo‘/src 


•WITH  AUTOPILOT  (SAS) 


*2*2$ +4 


HOVER  EXAMPLE  H 


WITH  AUTOPILOT 


NO  AUTOPILOT  (SAS 


TIME  SECONDS 


Figure  F-l.  Effect  of  Autopilot  Dynamics  Longitudinal 
Stability  1°  Longitudinal  Cyclic  Pulse. 


ROLL  ANGLE,  $ , 


TIME  r- SECONDS 


Figure F-2. Effect  of  Autopilot  Dynamic  Lateral  Stability 
Hover  1°  Lateral  Cyclic  Pulse  Example  H. 

f -\6 


i 


20 


AUTOPILOT 


NO  AUTOPILOT 


TIME  ~ SEC 


jF-3.  Effect  of  Autopilot  Dynomic  Lateral  Stability 
Velocity  s lOO'/sec  1°  Lateral  Cyclic  Pulse 
Example  I. 


Figure  f~L'  • Example  HI  Level  Flight  Trim  B vs  V Effect  of  Horizontal 
Tail  Coupling. 


CONTROL  LIMIT 

7777777 

C.G.  TRAVEL  (TOTAL) 
dit/dB,  « 0,  175" 
di*  /dBi  = I , 25.6" 


C.G.  TRAVEL  - INCHES 


Figure  F-5.  Example  IE  Estimated  C.G.  Travel  for  112 
Longitudinal  Cyclic  Control  Power  Trim 
Limits  Between  Hover  (Fwd  C.G.)  and 
High  Speed. 


a< 


LOAD  FACTOR  ~ 17 

Figure  F-7 . Example  Hr  Ml  Up  Velocity  100  kts 
Control  Position  vs  load  Factor. 


BANK  ANGLE~<£ 


Figure  F~8.  Example  HI  Control  Angle  vs  Bank  Angle 
Steady  Turn  V = 100  kts. 
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